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Abstract 

A  new  theory  for  extended  x-ray  absorption  fine  structure 
(EXAFS)  is  presented  which  treats  EXAFS  as  arising  from  the  scat¬ 
tering  of  the  ejected  photoelectron  by  atoms  surrounding  the  ab¬ 
sorbing  atom.  The  surrounding  atoms  are  treated  as  a  system  of 
weak  point  scatterers  which  allows  the  scattered  waves  to  be 
calculated  using  standard  scattering  theory.  These  scattered  waves 
are  included  in  the  wave  function  of  the  ejected  photoelectron  and 
this  is  used  to  calculate  the  dipole  transition  matrix  element 
from  the  x-ray  absorption,  which  gives  rise  to  oscillations  in 
the  absorption  coefficient  which  depends  on  the  positions  of  the 
surrounding  atoms.  The  final  result  includes  scattering  factor, 
range  and  temperature  terms  and  is  general  enough  to  handle  any 
symmetry  of  atoms  about  the  absorbing  atom  and  any  polarization 
state  of  the  x-rays.  The  result  for  unpolarized  x-rays  is  used 
to  calculate  EXAFS  for  Cu,  Fe  and  Ge ,  which  are  in  excellent  agree¬ 
ment  with  the  observed  spectra.  Furthermore,  the  results  of  the 
theory  can  be  used  to  show  that  when  EXAFS  is  Fourier  transformed 
it  yields  a  radial  structure  function  which  determines  the  dis¬ 
tribution  of  atoms  about  a  particular  kind  of  atom  in  the  absorb¬ 
ing  material.  A  method  for  determining  amorphous  structures  is 
developed  which  compares  the  amorphous  and  crystalline  radial 
structure  functions  and  uses  some  of  the  parameters  which  are 
determinable  from  the  crystalline  data.  It  is  shown  that  this 
technique  has  better  resolution  and  can  derive  more  information 
than  the  usual  radial  distribution  analysis.  Experimentally,  a 
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fully  automated  apparatus  is  described  which  takes  data  with  about 
.3%  statistical  counting  error  up  to  2000  eV  above  the  absorption 
edge.  Data  of  this  quality  is  necessary  in  order  to  obtain  reli¬ 
able  structural  information  from  EXAFS.  The  EXAFS  technique  was 
used  to  study  the  structure  of  amorphous  Ge,  Ge02,  and  GeSe.  In 
Ge,  the  data  was  consistent  with  a  continuous  random  network  of 
Ge  tetrahedra  where  the  order  decreases  very  rapidly  beyond  the 
second  neighbor  distance.  The  Ge02  structure  also  appeared  to  be 
a  tetrahedral  network  except  that  it  was  more  ordered  since  the 
adjacent  Ge04  tetrahedra  share  the  bridging  oxygen  atoms,  thus 
giving  a  correlation  between  adjacent  tetrahedra.  In  GeSe  the 
power  of  the  EXAFS  technique  to  determine  the  structure  about  each 
kind  of  atom  is  demonstrated.  The  crystalline  structure  is  a 
highly  distorted  NaCl  structure  with  a  3:3  near  neighbor  coor¬ 
dination.  The  amorphous  structure  shows  a  splitting  of  the  near 
neighbor  peak  indicating  a  tendency  for  both  Ge  and  Se  atoms  to 
satisfy  its  ideal  bonding  requirements.  The  results  of  this  work 
have  allowed  some  basic  conclusions  to  be  drawn  about  the  general 
nature  of  amorphous  structures. 


Chapter  1 
Introduction 


Extended  x-ray  absorption  fine  structure  (EXAPS)  refers 
to  the  oscillations  in  the  x-ray  absorption  coefficient  extend¬ 
ing  out  to  several  hundred  electron  volts  on  the  high  energy 
side  of  an  x-ray  photoelectron  edge.  This  phenomenon  has  been 
known  for  over  40  years  since  the  first  studies  of  Kronig1. 
However,  an  adequate  theory  which  quantitatively  explains  the 
observed  spectra  has  not  been  proposed  until  now.  Previous 
theories  have  been  reviewed  by  Paratt2  and  more  recently  by 
Azaroff^.  Pig.  1  shows  examples  of  EXAPS  for  the  Ge  K  edge  in 
amorphous  GeC^  and  for  both  the  Ge  and  Se  K  edges  in  amorphous 
GeSe.  Experimentally  In  I0/I  =  yx  is  plotted  as  a  function  of 
the  ejected  photoelectron  energy  where  I  and  IQ  are  the  x-ray 
intensity  with  the  absorber  in  and  out  of  the  beam,  y  is  the 
linear  attenuation  coefficient  and  x  is  the  thickness  of  the 
absorber. 

The  first  theory  of  EXAPS  was  proposed  by  Kronig-*-  who 
suggested  that  the  absorption  anomalies  (maxima  and  minima) 
were  caused  by  Bragg  reflection  of  the  ejected  photoelectrons 
from  the  planes  in  the  crystalline  lattice  causing  forbidden 
and  allowed  energies  for  the  electron  transition.  A  similar 
theory  was  proposed  by  Hayasi4  who  felt  that  the  Bragg  reflected 
electrons  led  to  oscillations  in  the  dipole  transition  matrix 
element  by  setting  up  standing  waves  (quasi-stationary  states) 
in  the  region  around  the  absorption  atom,  thus  increasing  the 
overlap  of  the  initial  and  final  electron  wave  functions.  Both 
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Examples  of  EXAFS  for  the  Ge  K  edpe  in  amorphous  0e02  and 
he  Ge  and  Ge  K  edpes  in  anorohous  GeSe.  The  last  two  curves 
have  been  displaced  downward  and  to  the  rinht  for  clarity. 


various  choices  of  scattering  potentials.  The  final  sro 
theory  was  proposed  by  Shiraiwa^,  et  al.  who  calculated 
scattering  from  individual  atoms  in  a  solid  including  multiple 
scattering  from  atoms,  and  included  a  range  term  for  the  photo¬ 
electron  due  mainly  to  inelastic  scattering.  They  concluded 
that  the  single  scattering  dominated  and  arrived  at  a  result 
which  was  similar  to  the  other  sro  theories. 

The  results  of  all  of  the  sro  theories  give  expressions 
for  the  fine  structure  which  are  equivalent  to 
X  v  sin  ^2kr  +  2n(k)^ 

j  r.  (1-1) 

J 

4“ 

where  r.  is  the  distance  from  the  absorbing  atom  to  the  j 
J 

coordination  shell,  N,-  is  the  number  of  atoms  at  r.,  k  is  the 
wave  vector  of  the  photoelectron  and  n(k)  is  the  phase  shift  of 
the  photoelectron  because  of  the  potential  of  the  absorbing 
atom.  Most  of  the  attempts  to  fit  the  data  have  involved  find¬ 
ing  the  peak  positions  of  the  dominant  shells,  and  the  agreement 

with  the  data  is  typically  within  5  -  10$.  Some  of  the 

5_q 

theories  have  added  scattering  factors-^  '  }  temperature 

O  Q 

corrections0,  and  range  terms ^  to  eq.  (1-1),  but  attempts  to 
calculate  entire  spectra  which  are  in  agreement  with  observed 
EXAFS  have  met  with  limited  success.  What  is  needed  to  calculate 
EXAFS  which  agree  with  experiment  is  a  simple  model  which 
retains  the  physics  of  the  problem  but  allows  straightforward 
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calculation  of  the  scattering  using  standard  pertubation  theory. 

A  photoelectron  range  term,  temperature  correction,  and  the 
energy  dependence  of  the  scattering  of  the  photoelectron  must 

all  be  included  in  any  quantitative  result. 

One  of  the  difficulties  in  comparing  theoretical  and 
experimental  results  is  that  early  experimental  data  was  not 
of  good  resolution.  The  first  experiments  placed  the  absorber 
in  the  x-ray  beam  before  the  analyzing  crystal  and  used  film 
to  record  the  spectra.  The  data  collected  in  this  way  is 
usually  of  poor  resolution  and  represents  only  the  variations 
in  the  intensity  through  the  absorber.  It  is  difficult  to 
divide  by  the  unattenuated  intensity  so  that  In  IQ/I  can  be 
plotted.  Another  limitation  of  this  technique  is  that  EXAFS 
were  recorded  over  only  a  few  hundred  eV  above  the  edge,  whereas 
EXAFS  can  be  observed  over  1000  eV  above  the  edge.  More 
recently  EXAFS  data  has  been  collected  using  scintillation 
counters  and  a  stepwise  recording  of  the  data.  Most  of  the 
experiments  have  been  manually  operated  so  that  the  data  record¬ 
ing  is  very  tedious  and  is  usually  limited  to  only  a  few  hundred 
eV  above  the  edge,  and  the  statistical  counting  error  at  each 
point  is  about  1 l.  In  order  to  obtain  EXAFS  of  sufficient 
quality  to  allow  comparison  with  theory  an  automated  experiment 
is  needed  which  can  gather  data  of  better  than  .3%  statistical 
counting,  error  over  an  interval  of  2000  eV  above  the  edge. 

Experiments  have  shown  EXAFS  to  be  sensitive  to  the 
atomic  structure  about  the  absorbing  atom.  However,  neither 
the  experimental  data  nor  previous  theories  have  been  sufficient 
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to  make  EXAPS  more  than  an  interesting  phenomenon.  In  this 
thesis  a  new  theory  and  improved  data  will  be  presented  which 
not  only  show  quantitative  comparison  between  theory  and 
experiment,  but  demonstrate  that  EXAPS  can  be  used  as  a  power 
ful  new  research  tool  to  study  the  structure  of  amorphous 


material . 


Chapter  2 

Point;  Scattering  Theory  of  X-ray  K  Absorption  Fine  Structure 


Any  successful  theory  will  have  to  explain  the  following 
experimentally  observed  conditions:  (1)  EXAFS  has  been  found 
in  crystalline  and  non-crystalline  solids,  liquids,  and  poly¬ 
atomic  gases  -  everything  except  monatomic  gases;  (2)  EXAFS  is 
sensitive  to  the  interatomic  distances  and  the  symmetry  of  the 
atoms  surrounding  the  absorbing  atom;  (3)  EXAFS  is  temperature 
dependent  and  becomes  sharper  at  lower  temperatures;  and  (^t) 
EXAFS  depends  strongly  on  the  symmetry  of  the  electron  that 
abscr-bs  the  x-ray.  For  instance,  the  and  spectra 

which  originate  from  p  electrons  are  identical  to  each  other 

but  different  from  the  and  K  EXAFS  which  are  identical  and 

10 

originate  from  s  electrons. 

In  this  chapter  a  simple,  flexible  theory  of  K  edge 
absorption  fine  structure  will  be  proposed  which  incorporates 
these  observations.  Although  the  physical  mechanism  responsible 
for  EXAFS  is  the  same  for  any  sro  theory  and  the  result  derived 
herein  have  elements  which  can  be  found  in  previous  sro 

O 

results,’  no  previous  alculation  incorporates  all  of  these 

necessary  elements.  There  are  also  differences  in  the  approxi¬ 
mations  that  are  made.  This  allows  calculation  of  spectra 
which  agree  with  the  observed  spectra  in  shape,  position,  and 
relative  size  of  the  peaks.  The  theory  does  not  attempt  to 
calculate  an  absolute  value  of  the  absorption  coefficient  but 
retains  only  those  terms  which  yield  an  oscillatory  behavior, 
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nor  does  it  calculate  the  structure  near  the  absorption  edge 
(i.e.,  within  50  eV  or  so)  which  is  affected  by  the  singularity 
produced  by  the  interaction  of  the  photoelectron  with  the 
ionized  atom,11  and  by  the  detailed  scattering  from  surrounding 

st oms * 

During  an  x-ray  absorption  event  conservation  of  energy 
demands  that  the  energy  of  the  ejected  photoelectron  ,E,  be 

E  =  hv  -  E  (2-1) 

K 

where  hv  is  the  energy  of  the  x-ray  photon,  is  the  K  shell 
binding  energy,  and  it  is  assumed  that  the  ion  core  is  left 
in  an  unexcited  singly  ionized  state.  If  Iq  is  the  intensity 
of  the  incident  x-ray  beam  then  the  intensity  of  the  beam 
transmitted  through  an  absorber  of  thickness  x  is  given  by  the 
familiar  equation 


where  y  is  the  total  linear  attenuation  coefficient.  Experi¬ 
mentally  I  and  IQ  are  measured  at  a  series  of  photon  energies 
and  In  (I0/I)  is  plotted  as  a  function  of  photoelectron  energy. 
Any  observed  fine  structure  is  attributed  to  oscillations  in 
the  linear  attenuation  coefficient.  This  coefficient  is  in 
turn  related  to  the  total  attenuation  cross  section  by 

y/p  =  o  (2-4) 

y  T 

where  a  is  the  number  density  of  atoms  in  the  absorber  and 
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°  =0+0  (2-5) 

T  A  sc 

where  o  is  the  photoelectric  absorption  cross  section  and n  i « 
A  sc 

the  total  scattering  cross  section.  a  will  be  ignored  since 

s  c 

it  does  not  contribute  to  the  oscillatory  behavior.  We  consider 
only  the  photoelectric  cross  section,  which  may  be  written  us¬ 
ing  the  dipole  approximation  for  radiative  transitions  12 

°A  =  1<1l2ahv  |e*£  |  2  N(E)  (2-6) 

where  a  is  the  fine  structure  constant,  N(E)  is  the  density  of 
final  states,  e  is  the  polarization  vector  of  the  x-ray,  and 
Mf  '*'S  d-^P°-*-e  matrix  element  given  by 


r 

•'•if 


r 

-v-2 


'^N;  ILj 


r 

^n 


n 


<!>  ( r 

f  Vl » 


£2 


(2-7) 

\VdT 


where  4>±  and  <p  are  the  total  many  body  wave  functions  of  che 
initial  and  final  atomic  states. 

To  evaluate  this  matrix  element  we  must  estimate  the  many 
electron  wave  functions  <j>^  and  which  are  given  approximately 
by  Slater  determinants  of  one-electron  wave  functions,^  and 
N-l  of  N  electron  states  which  are  the  same  for  the  initial  and 
final  states.  The  discrete  wave  functions  are  assumed  to  be 
normalized  to  unity  and  the  continuum  functions  normalized  per 
energy  interval.  (A  good  discussion  of  the  evaluation  of  this 


matrix  element  and  the  assumptions  relating  to  it  can  be  found 
12  13 

elsewhere.  5  )  Many  body  effects  are  now  included  as 


approximations  in  the  one  electron  wave  functions  which  are 
solutions  of  the  Hartree-Fock  equation.  However,  some  of  the 
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effects  of  the  hole  In  the  K  shell  for  energies  below  50  eV 
are  not  Included  In  the  wave  functions  such  as  those  that  produce 
singularities  at  the  edge11.  Since  shielding  can  distort  the 
wave  function,  particularly  at  energies  near  the  K  edge,  this 
theory  Is  not  expected  to  be  reliable  for  photoelectrons  of 


low  energy.  Under  these  assumptions  the  matrix  element  becomes 

m  -  /*{(?> dT  (2"8> 

Variations  In  the  absolute  square  of  this  matrix  element  are 
directly  related  to  the  x-ray  absorotlon  fine  structure  spectrum, 
in  the  case  of  K-shell  absorption  may  be  represented  by  the 
atomic  1-s  wave  function.  ♦,  Is  the  final  photoelectron  wave 
function  Including  scattering  and  It  Is  this  wave  function 
which  must  be  found. 


*f  Is  the  state  at  a  given  energy  which  Is  the  solution 
of  the  problem  of  a  solid  with  an  atom  at  the  origin  Ionized 
by  an  electron  missing  from  the  K-shell  while  the  surrounding 
atoms  are  neutral.  We  assume  a  muffin-tin-type  potential  - 
a  spherically  symmetric  potential  about  each  atom  with  a 
radius  which  Is  smaller  than  half  the  distance  between  atoms 
with  a  constant  potential  outside  in  the  rest  of  the  solid. 

The  ionized  absorbing  atom  will  have  a  different  potential  from 
that  of  the  neutral  atoms.  We  expect  the  potential  of  the 
ionized  atom  to  be  shielded  by  the  surrounding  conduction 
electrons  Justifying  Its  assumed  localized  nature.  The  ener¬ 
gies  of  interest  are  of  the  order  of  50  eV  or  so  greater  than 

These  energies  are  much  greater  than  any  of 


the  Fermi  energy. 


11 


the  expected  matrix  elements  of  the  pseudo-potentials  or  T- 
matrices  that  determine  the  scattering  from  the  neutral  atoms. 

For  this  reason  the  scattering  should  be  weak  and  contributions 
which  correspond  to  multiple  scatterings  from  the  neutral  atoms 
can  be  neglected.  Therefore,  each  neutral  atom  can  be  replaced 
by  a  delta  function  potential  chosen  to  give  the  same  back- 
scattered  wave  toward  the  ionized  atom  as  would  the  actual 
neutral  atom. 

The  ionized  atom  is  more  difficult  to  treat.  It  cannot  be 
replaced  by  an  effective  potential  since  overlap  of  the  initial 
and  final  state  wave  functions  must  be  determined  near  the 
nucleus.  Hence,  the  full  problem  must  be  solved  with  the 
complete  potential.  The  variations  in  the  x-ray  absorption 
spectrum  are  produced  by  the  interference  near  the  nucleus  be¬ 
tween  the  outgoing  wave  from  the  absorbing  atom  and  the  waves 
backscattered  from  the  surrounding  neutral  atoms.  When  these 
waves  add  constructively  there  is  an  absorption  maximum  and  a 
minimum  when  they  add  destructively.  The  problem  of  solving  the 
final  state  wave  function  inside  the  ionized  atom  will  be  de¬ 
layed  until  the  form  of  the  wave  function  Just  outside  of  the 
ionized  potential  is  considered  (where  the  potential  is  constant). 
The  wave  function  will  have  two  parts,  an  outgoing  wave  corres¬ 
ponding  to  an  electron  excited  from  the  K-shell  and  a  part 
consisting  of  backscattered  waves  from  the  surrounding  neutral 
atoms.  The  outgoing  wave  can  be  separated  into  radial  and 
angular  parts.  In  general,  the  radial  part  will  have  the  form 
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R  ■  Nh^(kr)  (2-9) 

where 

h'(kr)  =  exp  (2in  )h*(kr)  (2-10) 

and  N  is  si  normalization  constant  assumed  to  be  real  ,  h^  is 
the  outgoing  Hankel  function  and  may  be  visualized  as  an  out¬ 
going  spherical  wave  relative  to  the  ionized  atom.  The  phase 
shift  is  introduced  by  the  presence  of  the  ionized  atom  at 
the  origin  and  would  be  zero  if  there  were  no  such  atom  present. 

For  K-shell  absorption  the  initial  electron  state  is  an  s 
state  and  the  dipole  selection  rule  [Afc  =  ±1  for  the  matrix 
element  (eq.  2—8)]  allows  only  transitions  to  states  with  p 
symmetry  ( 8.  ■  1).  The  total  wave  function  will  be  the  product 

of  h'(kr)  and  a  linear  combination  of  the  three  spherical 
1 

harmonics  for  fc  =  1  which  represent  the  x,  y,  and  z  components 
of  the  wave  function.  These  components  will  be  considered  one 
at  a  time  starting  with  the  z  component.  The  spherical  harmonic 
in  t£he  z  direction  is  Y°  =*  cosQ,  so  that  the  wave  function  of 
the  outgoing  part  of  the  photoelectron  wave  function  is 

i|»(r)„  *  Nh,’  (kr)  COS0  (2-11) 

Z  X 

To  obtain  the  total  wave  function  we  must  add  the 
contribution  from  the  scattered  waves.  This  contribution  is 
composed  of  the  scattering  from  all  of  the  surrounding  atoms. 

As  discussed  previously,  we  are  interested  in  the  energy  range 
where  scattering  from  a  single  atom  is  small.  We  will  there¬ 
fore  neglect  the  multiple  scattering  contribution,  i.e., 
scattering  from  more  than  one  atom.  However,  we  will  include 
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multiple  scattering  from  the  same  atom  which  converts  the 
potential  matrix  to  the  T-matrix.  It  is  only  the  T-matrix 
which  is  small  from  an  isolated  atom.  The  scattered  wave  can 
now  be  calculated  by^ 


sc 


5^(E-H  -A7]T(r  )  **  (2-12) 

*"T  o  '''if 


where  (E-H  -A)”1  is  the  Green's1  function  of  the  average  medium 
o 

for  energy  E,  through  which  the  electron  propagates.  T^)  is 
the  T-matrix  of  the  neutral  atom  at  ,  and  is  the  incident 
wave  at  point  r^  produced  by  ail  other  neutral  atoms  excluding 
the  atom  at  r^  and  is  given  by 

*r  ’  *  +  E 

mj*i 

and  A  is  an  operator  with  only  diagonal  elements  equal  to  nTkk 

defined  by  eqs.  (2-16)  and  (2-17).  In  our  approximation  of 

small  T-matrix  ^  is  independent  of  the  atom  at  r^  and  is  equal 

to  ii*„.  We  replace  T(r  )  by 
f  'vm 

(2-14) 


(E-H  -A)"1  T(r  )  (2-13) 

o  ^m  f 


TtCm>  *  Vo  5(^m> 


with  V  chosen  to  produce  the  same  backscattered  wave  as 
o 

produced  by  T(^ft).  The  equation  we  need  to  solve  is,  substitut¬ 


ing 


and  using  eq.  (2-l4) 


♦so  ■  vo  E  (E-vA)"1 


(2-15) 


It  Is  important  to  note  that  is  significantly  different  from 
iji.  As  ij<f  moves  further  from  the  origin  more  and  more  of  the 
scattered  waves  from  surrounding  atoms  are  added  to  it  and  it 
deviates  more  from  <J*.  Sufficiently  far  from  the  origin  <jif 
should  be  modified  to  the  form  of  an  eigenstate  and  produce  no 
more  scattered  waves.  In  addition,  electron— electron,  and 


electron-phonon  scatterings  make  ^  lose  coherence  with  *  . 

Thus,  only  the  atoms  within  a  certain  distance  of  the  orieln 
contribute  to  the  backscatterlng.  An  estimate  of  this  scattering 
length  will  be  given  later.  This  effect  is  approximated  by 
assuming  that  *  decays  exponentially  away  from  the  origin. 

The  scattering  has  two  other  effects.  The  scattering  from 
the  periodic  array  of  neutral  atoms  produces  energy  gaps  cor¬ 
responding  to  energy  regions  where  no  propagating  states  can 
occur.  In  these  regions  the  E  vs .  k  spectrum  and  the  density  of 
states  are  distorted  from  a  simple  free  electron  behavior. 
However,  because  the  scattering  Is  weak,  the  widths  of  the  gaps 
are  small,  and,  since  the  experimental  data  averages  over  several 
eV  which  is  an  energy  Interval  large  compared  to  the  gaps ,  the 
gap  effects  should  not  be  observable.  For  that  reason  we  neglect 
the  gap  effects  on  the  density  of  states,  and  assume  that  11(E) 
varies  as  in  a  free  electron  gas,  proportional  to  k.  The 
scattering  from  the  periodic  array  also  changes  the  average  k- 


(2-16) 


value  of  41  to  k'.  The  equation  determining  this  is 
f 

2  2 

E  -  - n  T  =  0 

2m  k'k' 

where  n  is  the  number  of  neutral  atoms  per  unit  volume  and  T 
is  the  forward  scattering  amplitude  Riven  by 


k'k' 


k'k' 


fexp(-11&' exp(ij£*  •£)  dx  (2-17) 


Equation  (2-15)  can  be  written  as 


V  2m 

(r)  - - - — 

sc 


(2-18) 

V  r  exp(ik’  Ie-e'  I  )  6(r_r.)  (r«)  dx 

L  J  -  ^  n  f  * 


j  r-r  ’ 


where  the  spatial  representation  of  the  Green  s  functi 
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been  used  and 


4V 

2m 


As  discussed  above  Is  approximated  by 

d/-(r)  ■  Nhi  (k'r)  cose 
i  z  1 

where  k'  Is  determined  from  eq.(2-26)and 
scattering  length  mentioned  previously. 


and(2-20)we  obtain 

V  2mN 

’ll  (r)  =  -  _2 - 


sc 


2nfi 


sr  expdk'l^-^l  ) 
lr-ril2 


(2-19) 

exp(-yr)  (2-20) 

y  Is  a  measure  of  the 
Combining  eqs.(2-l8) 

h'(k'r,)(z-z  )  exp(-Y r^) 
(2-21) 


where  the  fact  that  the  cosine  of  the  angle  between  r  and  the 
z-axls  is 


(z-z  ) 
cos  e  ■  _ L. 

1  lc-?1 1 


(2-22) 


arid  z^  is  the  z-component  of  r^  has  been  used.  The  total  wave 
function  outside  the  Ionized  atom  will  be  a  combination  of  eqs. 
( 2-11) and(  2-211  giving 


(2-23) 


*f(r) 


/  x  v  mN  exp(ik'|  r-r  |)  .....  »  -  yr .  .  . 

Nh'.(k'r)  cose  _  >  _  1  h' (k'r  )  e  1  (z-z  ) 

1  _  L*  -  11  1 

I r-r  |2 


2rrft 


This  expression  for  is  expected  to  be  a  reasonable 
approximation  just  outside  the  Ionized  atom.  It  remains  to 
extrapolate  this  expression  inside  the  ionized  atom.  This  is 
done  by  using  the  WKB  approximation  with  the  Langer  correction  for 
the  phase  shift.  This  phase  shift,  introduced  by  the  presence  of 
the  ionized  atom,  is  given  by1^ 
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(2-24) 


ni = 


r 


2mV(r) 
2 


•ft 


2  25\1/2  f°/  ,2  2.25\1/2 

-f)  dr-  /(k  'V 


dr 


where  r  and  r  are  the  zeros  of  their  respective  integrands,  V(r) 
is  the  atomic  potential,  and  ro  is  a  point  on  the  constant  cart  of 
the  rauffin-tin-potential  (assumed  to  be  zero).  In  determining  * 
near  r  =  0  this  phase  shift  must  be  subtracted  from  the  first  term 

and  added  to  the  second  term  on  the  right  side  of  eq.  (2-23). 

For  the  medium  and  heavy  elements  the  1-s  function  Is  snarply 
peaked  near  the  oriein  and  the  o-wave  can  be  approximated  by  a 
straight  line  passing  through  the  origin,  thus  the  amount  of  o- 
character  in  the  wave  should  be  closely  approximated  by  the 
magnitude  of  the  slope  of  this  line.  The  matrix  element  eq  (2-8) 
will  therefore  be  proportional  to 


M  6  — - - 

^  'sry 


N 


3  h+(0)cos  9 
1  _ o 

az 


exo(2in)V  mn 
‘ _ 1  o 

2  -rr+i2 


(2-25) 


x 


a__ 

az 


I 


exp(ik’ j  r_ri !  ) 


|  r-ri l 


2 


h+(k’r  )  exp(-yr  )  (z-z  ) 
1  i  11 
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(2-26) 


M  6  A  +  C  exp(2in  )V  A 
%  1  o 


exp(ik'r  ) 

- = -  h  (k*r  )  exp(-yr  ) 

r  2  11  1 
1  i 


2z 


2 


ik'z 


-  1  - 


where  A  is  a  constant  proportional  to 


ah  (0) 

1 

az 


The  x-ray  absorption  cross  section  per  unit  energy  is  pro¬ 
portional  to  I ( e • M ) | 2  k',  where  the  k*  comes  from  the  density  of 
final  states  eq.  (2-6).  Only  the  cross  product  in  |M|2  contributes 
to  the  oscillatory  fine  structure  spectrum.  Therefore,  the  fine 
structure  spectrum  is  proportional  to 


X  =  Re 


exp(ik'r  ) 

V  1  + 

e  2  exp(2in  )V  k'  >  - 2 -  b  (k'r  )  exp(-yr  ) 

z  1  o  11  i 


i 

*  <7.  ^  \  1 


(2-27) 


where  Re  [  ]  means  the  real  part  of  the  bracketed  term.  Using 
the  relationship 


h+(kr)  =  -  exp(ikr) 
1 


(t*r)2  +  to) 


(2-28) 


in  eq.  (2-27)  we  can  write 
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(2-29) 


where 


a 


1 


sin[( 2k '  )r  +2n  )] 
_  1  1 


k'r 

1 


cos[(2k,)r1  +  2n1] 


(2-30) 


and 


8 

1 


cos[(2k')r  +  2n  )] 

1  1 

- — - —  + 

k'r 

1 


sln[(2k' )r^ 


+  2n  ] 


(2-31) 


Equation  (2-29)  may  be  simplified  by  comparing  the  relative 
magnitudes  of  the  separate  terms  and  noting  that  experimentally  k' 
varies  from  about  2  to  20  A  ,  while  r  varies  from  about  1.5ft  to 
10ft.  Using  these  numbers  and  simplifying;  the  notation  by  dropping 
the  crime  on  k  in  eq.  (2-29)  reduces  to 


2  yexp(-Yr  ) 

x  %  -Vok  ez  Z-— p j -  Zj2  sin[2kr1  +  2n(k)]  (2-32) 

1  ri 

It  now  remains  to  calculate  the  fine  structure  contribution 
of  the  x  and  y  components,  however  the  task  is  considerably  simpli¬ 
fied  by  noting  that  both  the  derivation  and  the  final  expression 
(eq.  2-29)  are  invariant  under  rotation  of  the  coordinate  axes  so 
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that  the  final  answer  may  be  written  directly 


exp(-yr  ) 


x(k)  •  -V  k  4,  5" 

°  i  ri 

/  2  2  2  2  2  2  \ 
(ez  zi  +  ey  yi  +  ex  xi  ) 


sin^2kr1  +  2n(k)>) 


(2-33) 


This  expression  may  now  be  evaluated  for  various  crystal  sym- 
metries  and  polarization  states  of  the  incoming  x-rays  however  the 

2  C.  £  _  ^  /g 

case  of  interest  is  unpolarized  x-rays  so  thau  e^  e^  z 
which  gives  the  final  result 


x(k) 


p  p 

kVo  Yje  J  “2o-  k 

= .  —  z 


-^O  K  .  . 

■ 2 -  e  ^  sin(2kr  +  2n  )  (2-3*0 

j  rj  J  1 


where  the  sum  is  now  over  coordination  shells,  is  the  number  of 
atoms  in  the  jth  shell,  the  term  exp(-2^V)  U  a  temperature  ^ 

term  to  be  discussed  later,  and  x*  +  y*  *  **  *  ^  .  Eq.  (2-33) 
may  also  be  evaluated  for  materials  with  cubic  symmetry  for  any 

polarization  of  the  x-rays.  Since  x^  -  ri  S3 

cubic  symmetry,  eq.  (2-33)  reduces  to  eq.  (2-3-0  for  this  case  also. 

*  _  n  ll  \ 


EXAFS  spectra  calculates  with  this  simplified  expression  eq.(  3 
have  changes  in  peak  positions  of  less  than  1  eV  and  changes  in 
in  amplitude  of  less  than  2%  from  the  spectra  calculated  with  the 
complete  equation  eq.(2-29).  The  importance  of  the  separate 
factors  in  this  equation  for  calculate  EXAFS  spectra  will  now  be 

discussed. 
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Note  that  the  photoelectron  wave  vector  k  is  not  related  to 
E  by  the  simple  free  electron  relation  shown  in  eq.  (2-19),  but 
instead  by  eq.  (2-16)  where 


nT  -  (W  -  W) 
kk  o 


(2-35) 


and  W  is  the  work  function  and  Wq  is  the  average  potential  energy 
of  an  electron  at  rest  between  the  atoms  in  a  crystal.17  since 
the  experimental  data  is  taken  as  a  function  of  energy  while  the 
theoretical  expression  is  a  function  of  k  it  is  important  that  a 
reliable  method  of  calculating  or  determining  W  be  found. 
Kozlenkov  has  calculated  Wq  for  a  number  of  elements  by  averaging 
;he  atom's  potentials  over  a  unit  cell  of  the  crystal.  Mott17 
laa  devel°Ped  a  method  in  which  he  arbitrarily  chooses  W  and 
references  the  experimental  spectrum  to  this  new  energy  scale.  He 
uses  eq.  (2-16)  to  find  the  value  of  the  extrema,  and  finds 
phase  shifts  of  these  extrema,  and  compares  them  to  the  WKB 
hase  shifts  he  has  calculated.  This  process  is  repeated  with  a 
®W  choice  of  WQ  until  the  best  fit  with  the  calculated  phase 
bift  is  obtained.  For  general  application  this  method  has  the 
isadvantage  that  the  WKB  phase  shifts  may  not  be  the  best  phase 
lifts  to  use.  This  depends  on  the  choice  of  potentials  and  in 
my  materials  suitable  solid  state  potentials  are  not  known  and 
;ber  less  satisfactory  potentials  must  be  used.  Also  the  EXAFS 
iat  Mott  measured  originated  primarily  from  the  first  coordination 
'ell.  In  many  materials  the  observed  EXAFS  has  major  contri¬ 
tions  from  several  shells  and  it  is  difficult  to  sort  out  which 
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features  belong  to  which  shells  over  the  entire  range  of  the 
observed  spectrum.  On  the  other  hand,  calculations  like  those  of 
Kozlenkov  are  tedious  and  are  appropriate  only  for  crystalline 
materials  where  the  structure  is  known.  This  problem  is  related 
to  the  general  uncertainty  in  the  solid  state  potentials  for  the 
wide  variety  of  materials  which  can  be  studied  by  EXAFS,  i.e. 
elemental  substances,  alloys,  and  amorphous  metals  and  semicon¬ 
ductors.  If  the  potential  is  unknown  then  both  the  internal 
potential  and  the  phase  shift  are  uncertain  and  some  method  which 
finds  the  quantities  self-consistently  to  obtain  a  "best  fit  of 
eq.  (2-3*0  to  the  observed  spectrum  must  be  found.  For  simple 

metals  such  as  Cu  and  a-Fe  reasonable  potentials  exist  so  that 

reliable  phase  shifts  can  be  calculated.  The  inner  potentials 

are  usually  reliable  to  within  a  few  eV  and  a  good  fit  to  the 

observed  EXAFS  can  be  obtained  by  varying  the  value  of  the  inner 

potential.  This  has  the  effect  of  shifting  the  energy  scale  by 

the  amount  of  the  variation  and  expanding  or  contracting  the 

spectrum  slightly  by  changing  the  periodicity  of  the  sine  waves 

in  eq.  (2-3*0.  For  more  complicated  materials  a  variational 

method  similar  to  Mott's  using  just  the  nearest  neighbor  distance 

in  eq. (2-3*1)  can  be  used  provided  that  reliable  phase  shifts  and 

near  neighbor  distances  are  known.  Although  finding  reliable 

values  for  W  -  W  is  essentially  an  unsolved  problem  for  complex 
o 

materials,  the  shapes  of  spectra  calculated  from  reasonable 
estimates  of  this  quantity  are  not  sensitive  to  variations  in 
(which  only  shifts  the  spectra  and  varies  the  positions  of  the 
extrema).  This  does  not  alter  the  basic  correctness  of  our  result. 
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The  factor  VQ  in  eq.  ( 2—3^ )  is  chosen  to  give  the  correct 

backscattering  of  the  photoelectron.  In  backscattering  the 

electron  changes  its  wave  vector  by  2k  and  thus  V  <6  V(2k)  the 

o 

2kth  Fourier  component  of  the  scattering  potential.  This,  in 
turn,  is  related  to  the  electron  scattering  factor,  f,  and  hence 
forth  the  identification 

V  =  f ( 2k)  (2-36) 

o 


will  be  made,  which  is  expected  to  be  a  reasonable  approximation 
of  the  k  dependence  of  the  scattering.  A  qualitative  understanding 
of  the  z  dependence  of  the  backscattering  can  be  obtained  from  the 
following  argument.  Physically  the  backscattering  must  be  due 
to  scattering  from  the  ion  core  in  order  to  produce  the  needed 
large  momentum  transfers.  The  amount  of  penetration  of  the  photo¬ 
electron  into  the  ion  core  may  be  estimated  by  noting  that  back- 
scattering  occurs  from  the  charge  within  a  region  on  the  order 

A 

of  twice  the  wave  length  of  the  photoelectron  (i.e.  d a around 
the  scattering  atom.  The  ionic  charge  that  this  represents  can 
be  estimated  from 


Z  =  Z  -  n 
eff  eff 


(2-37) 


where  Z  is  the  atomic  number  of  the  atom  and  n  __  is  the  number 

eff 

of  electrons  within  the  radius  d  around  the  scattering  atom.  This 

may  be  estimated  from  the  electronic  structure  of  the  particular 

atoms  involved  and  the  radii  of  the  electronic  shells  given  by 
18 

Slater  .  The  result  of  these  estimates  for  0,  Ge,  and  Se  for 


some  typical  values  of  k  are  given  in  Table  1. 
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TABLE  1 


0 


Ge 


Se 


Data  for  Back-Scattering  in  0,  Ge  and  Se 


kdf1) 

E(eV) 

d(A) 

Z 

eff 

fo0O 

k  f(k) 
zeff 

5 

95 

.63 

1 

3.8 

95 

10 

380 

.31 

5 

2.4 

48 

15 

1520 

.16 

6 

1.7 

63 

k(A_1) 

E(eV) 

d(ft) 

Z 

eff 

f  (k) 

0 

2 

k  f (k) 

Zeff 

5 

95 

.63 

4 

9.7 

60 

10 

380 

.31 

13 

8.0 

62 

15 

1520 

.16 
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6.4 

65 

k(jf  ) 

E(eV) 

d(A) 

Z 

eff 

f  (k) 

0 

z 

eff 

5 

95 

.63 

6 

10.3 

43 

10 

380 

.31 

15 

8.3 

55 

15 

1520 

.16 

24 

6.7 

63 
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E  is  the  photoelectron  energy  for  the  listed  values  of  k 

calculated  from  eq.  (2-19).  A  comparison  of  the  Z  and  the 

eff 

19 

scattering  factors  from  Doyle  and  Turner  are  also  given,  and 

since  the  scattering  from  a  point  charge  of  Z  is  expected  to  be 

eff 

Z  e 

f(k)  =  4,  -  (2-38) 

2 

then  k  f(k)/Z  should  be  approximately  constant.  It  can  be  seen 
eff 

2 

from  Table  1  that  the  values  of  k  f(k)/Z  particularly  at  higher 

energies  are  approximately  constant  and  equal  to  about  the  same 

value  for  all  atoms.  This  implies  that  f(k)  scales  approximately 

as  Z  as  required  by  eq .  (2-38).  From  Table  1  the  values  of 
eff 
Z 

eff  are  roughly  proportional  to  Z.  This  allows  an  average  f  to  be 

defined,  as  a  rough  approximation  for  a  multicomponent  material, 

which  is  proportional  to  the  average  Z  for  the  material.  Although 

it  is  much  better  to  use  the  tabulated  f(k)'s  for  each  component 

the  present  quality  of  the  data  is  such  that  the  general  behavior 

of  the  backscattering  can  be  described  by  an  average  f(k). 

The  product  kf(2k)  forms  the  envelope  of  the  fine  structure 

spectrum.  In  Fig.  2  the  product  kf(2k)  is  compared  for  a  free 

electron  scattering  factor,  a  screened  coulomb  scattering  factor, 

19 

and  the  self-consistent  scattering  factors  of  Doyle  and  Turner. 

All  three  curves  give  a  similar  behavior  at  high  energies.  At  low 
energies  the  free  electron  factor  diverges  while  the  other  two 
curves  give  a  peak  about  k  =  4  and  fall  to  zero  as  k  •*  0  which  is 
better  agreement  with  the  observed  EXA^S.  Since  the  theory  isn’t 


Figure  2 

Comparison  of  the  factor  kf(k)  using  (a)  the  free  electron 
(b)  the  screened  coulomb,  and 
(c)  the  Doyle  and  Turner  scattering  factors. 
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expected  to  be  valid  below  about  50  eV,  either  scattering  factor 
is  acceptable.  The  scattering  factors  of  Doyle  and  Turner  are 
especially  convenient  since  they  are  presented  in  terms  of  an 
expansion  in  a  ^  term  exponential  series  which  lends  itself  to 
direct  computer  coding. 

The  range  term  exp(-yr)  in  eq.  (2-3*0  is  similar  in  form 

g 

to  that  used  by  Shiraiwa  .  However,  he  attributed  the  range  to 
inelastic  scattering  which  is  questionable  since  estimates  of  this 
scattering  length  are  on  the  order  of  several  lattices  spacings 
for  the  photoelectron  energies  of  interest.  This  is  an  apparent 
contradiction  to  experiment  which  shows  only  significant  contri¬ 
butions  from  about  two  lattice  spacings  (^5  coordination  shells). 
Instead  we  propose  that  the  dominant  operative  as  the  photo¬ 
electron  propagates  farther  and  farther  into  the  lattice  is  that 
the  scattered  waves  from  the  surrounding  atoms  are  added  to  the 
wave  function  and  it  becomes  more  like  an  eigenstate  of  the  system, 
and  therefore,  produces  no  more  coherently  scattered  waves,  A 
rough  estimate  of  this  distance  can  be  made  as  follows :  The 
distance,  d,  that  the  photoelectron  travels  before  it  is  scattered 
into  an  eigenstate  of  the  system  can  be  estimated  from 


m 

where  t  is  the  time  it  takes  to  be  scattered  into  an  eigenstate 
and  is  not  related  to  the  usual  lifetime  due  to  inelastic  scattering. 

t  can  be  estimated  by  considering  the  expansion  of  the 
electron  wave  function  of  the  solid 


<Kr,t) 


a  ip 
n  n 


n 


(2-40) 
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where  the  sum  is  over  the  states  of  the  system  which  differ  by 
reciprocal  lattice  vectors  and  is  subject  to  the  condition 


4,(r,0)  =  exp(ik*r)  (2-Mi; 

Generally  in  a  solid  the  two  dominant  states  in  the  sum 

will  be  k  and  k  +  p,  where  p^  is  the  shortest  reciprocal  lattice 

vector.  Considering  only  these  two  states  the  phase  term  in 

en  (2-40)  is  exp[i(E  -E  )t/ft].  This  may  be  rewritten  using 

k 

the  definition  for  the  free  electron  energy  eq.  (2-19)  and  the 

fact  that  k  »  r,  for  the  electron  energies  of  interest.  Then 

min 


r1(\+Ri  -  V‘‘ 

a  exp 

Iftkg  t" 

exp 

.  m 

When  this  argument  is  non-zero  the  plane  wave  term  will  be 
important  and  the  other  terms  will  partially  cancel  each  other; 
however,  when  the  argument  of  the  exponential  is  n,  all  terms 
add  directly  and  change  the  wave  function  to  an  eigenstate  of  the 
system.  This  gives  an  estimate  of  the  time,!. 


T 


mr 

tfkg 

1 


(2-43) 


In  the  case  of  copper  which  has  a  b.c.c.  reciprocal  lattice 


P, 

1 


(2-44) 


Therefore,  from  eq.  (2-39)  and  (2-W) 


d  =»  a//T 


(2-^5) 


or  d  -  2.1*  In  copper.  This  mechanism  dominates  all  other 
scattering  mechanisms  in  determining  the  range  over  which  the 

Photoelectron  backscattering  contributes  to  EXAPS.  The  result  is 
also  independent  of  k  for  k  „  g^  Thls  arRunent  ldentmes  ^ 

range  parameter  Y  as  the  reciprocal  of  the  scattering  length  for 

scattering  of  the  photoelectron  into  an  eigenstate.  Its  accurate 

determined  by  a  fit  to  experimental  data  for  EXAPS  and 

it  typically  about  v  -  .3.  In  Chapter  A,  Y  will  be  determined 

directly  from  the  experimental  data  for  crystalline  germanium. 

The  term  exp^2*?)  ls  a  Beneral  Debye  form  fop  the 

temperature  correction  that  must  be  applied  because  of  the  thermal 
motion  of  the  lattice. 

The  origin  of  this  term  may  be  seen  by  the  following  simple 
derivation.  If  the  thermal  motion  of  all  of  the  atoms  in  the 

solid  is  taken  into  account  then  the  theoretical  expression  for 
EXAPS  should  be  written 


x(t)  =  -kf<k>2  sin(2kr,(t)  ♦  2„, 

r*  c  J- 


(2-1*6) 


where  r^t)  is  the  distance  between  the  absorbing  atom  and  the 
1th  atom,  including  the  thermal  motion  of  both  atoms,  and  the  sum 
is  now  over  Individual  atoms.  It  Is  the  time  average  of  eq.  (2-« 
that  Is  needed  for  EXAPS.  If  one  typical  sine  term  is  considered 
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from  the  sum  then  the  time  averaged  expression  can  be  simplified 
by  rewriting 


r^(t)  =  r^°  +  a^(t) 


(2-47) 


where  r^0  is  the  equilibrium  distance  between  the  absorbing  atom 
and  the  atom  and  o^(t)  is  the  relative  displacement  of  the 
atoms  because  of  the  thermal  motion  then  we  may  write 


sin(21cr^lt)  +  2nJ  =  sin(2kr^°  +  2n)  cos  ^ko^Hf)  jj- 


(2-48) 


cos(2kr^°  +  2n)  sin  ^2ko^  ( t )^ 

where  the  bars  indicate  time  averaging.  Since  the  disnlacements 
are  random  in  sign  and  magnitude  the  time  average  of  sin|2ko^(t)) 
will  be  zero  so  that 


sin(2kr.(t)  +  2n)  =  sin(2kr°  +  2n)  cos(2ko  (t) 
i  i  i 


(2-49) 


the  cosine  term  may  be  expanded  as 

2 

cos  ^2ko^(t )|-Ci-  1  -  (-2-k— •  a^(t)  + 


(2-50) 


but  this  can  be  related  to 

2 

-x 

e  .ci 


2 

1  -  x  +  . . . 


(2-51) 
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where 


2  2  2,  % 

x  =  2k  (t) 


(2-52) 


so  that  the  expression  for  EXAFS  becomes 


N  e 

x (t )  %  -kfoo  y  2 


-Yr  2  2 

J  -2k  o4(t)  ,  „  . 

“  e  .1  sin(2kr  +  2n) 

j 


(2-52) 


where  the  sum  is  now  over  coordination  shells  and  it  is  assumed 
a  2(t)  is  the  same  for  all  atoms  in  the  coordination  shell. 

J 

Although  this  derivation  was  simplified  to  show  the  origin  of  the 

temperature  term  in  EXAFS;  it  can  be  made  more  rigorous  and  shown 

21,22 

to  be  valid  over  a  large  temperature  range  .  The  shell  to 
shell  variation  in  o  can  be  seen  by  relating  a  to  the  individual 
motion  of  the  absorbing  and  scattering  atom,  as  follows 


o  2(t)  =  u°  (t)  +  u  (t)2  -  2u°(t)  u  (t) 

J  J  J 


(2-5^) 


where  u  (t)  and  u  (t)  are  the  displacements  of  the  absorbing  and 

J  J 

scattering  atom  along  the  equilibrium  position  vector  and  this 

2 

expression  is  valid  to  order  /  ' 


(r 


If  u  ■  u  as  expected  in  a  monatomic  material  then  eq.  (2-54) 
j  J  p  2 

shows  that  a  »  2u  for  two  atoms  vibrating  completely  independ- 

.  j  j  — - - -  2 

ently  of  each  other  (i.e.  u°(t)  u  (t)  =  0)  while  o  =0  for  atoms 

J  .1  _  J _  2 

q  a  1  1  C. 

whose  motion  is  completely  coupled.  (i.e.  (u°(t)  u(t)  =  2u  ). 

J  j 
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An  exact  theory  for  the  temperature  dependence  of  any  material 
depends  on  calculating  u?(t)  u  (t)  for  each  shell  which  in  turn 

«]  iJ 

depends  on  the  structure  and  force  constants  in  a  material.  This 
does,  emphasize  that  the  temperature  correction  will  he  different 
Tor  each  shell  in  a  material  chan/rinp;  from  a  more  couoled  behavior 


•*M  “ije  first  shell  to  the  independent  value  in  the  higher  shells. 

The  general  behavior  of  a  for  different  shells  nay  be  seen 

by  writing  the  u's  in  terns  of  the  normal  modes  of  the  lattice,  thus 

J 


tTjTn; 


ft  (<)  ^A<)  exp(iK*r(m)) 

o  ^  ni  %  -\,  ^ 


wuere  n  is  the  index  for  the  three  nhonon  polarization  modes,  < 
is  tne  phonon  wave  vector;  afl  is  the  amplitude  of  the  normal  mode 
\  is  the  unit  disnlacement  vector  for  the  mode,  m  is  the  mass  of 
u!,e  ai-ons  (monatomic  material)  and  hT  is  the  number  of  atoms  and 
the  time  average  has  been  taken.  Usinn  en.  (2-55)  we  may  write 


u1(in)u1(n)  =  ^  (*)  *%,(«')  {  exn(i(K  •r(m)-^'  ^(n)  )j 


Reproduced  from 
_bes>  available  copy. 


(2-56) 


this  expression  is  for  one  pair  of  lattice  points.  It  must  be 
a ve raped  over  pairs  of  lattice  points  separated  by  r(m)-r(n)  so  that 


X  exp(iK»r  )  exo( 1 («-<’)• r(n) ) 
^  ^mn  -v  %  % 


(2-57) 
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where  r  =  r(m)-r(n),  but 
'v-mn  %  ^ 

2  exp(i (k-k ' ) • r (n) )  =  NA  (2-53) 

^  ■'V  .  -v  1C  j  1C  * 

n 


so  that  eq.  (2-57)  becomes 


Ui (m)ui (n) 


2  !%<J£>  I2 


exp(iic 

'V/ 


r  ) 
^mn 


(2-59) 


where  we've  used  the  fact  that  |  is  a  unit  vector  and  that  if 
«=<'  then  o=a '  since  the  modes  are  indenendent .  Eq.  (2-59)  and 
(2-5^)  may  now  be  combined  to  qive 


y  |a  (<)|2  (l-exp(iic  .r  )) 
/L  a  *  *  ^nn 

tc  ,  o 


(2-60) 


At  this  point  we  start  approximations.  Let  us  assume  the  Debye 

approximation  where  the  energy  of  a  mode  depends  only  on  the 

magnitude  of  k  for  a  Riven  o.  It  then  follows  that  a  (*)  denends 

a  ^ 

only  on  the  magnitude  of  c.  By  changing  the  sum  over  k  to  an 
integral  and  doinm  the  integral  over  angles 


on  '  mn  ' 


' 2  die 


(2-61) 


o  o 

where  <  is  the  maximum  phonon  wave  vector;  and  V  is  the  volume  of 
the  solid.  Tiie  amplitude  of  the  normal  waves  can  be  expressed  in 


terms  of  the  energy  of  the  wave  usinp;  the  usual  expression  for  the 
harmonic  oscillator 


|a  U)j2 


2fi 

mNa)  ( ic ) 


(n  +  i) 


(2-62) 
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where  n  is  the  Bose-Einstein  distribution  r.iven  by 


n  =  (exp^-ll)-!)"1  (2-63) 

1  k  T  ' 

B 


thus  eq.  (2-6l)  becomes 


1  _V_  2h 

N  2* 2  n 


fmC(e*U/ V  -I)"1 


+ 


sin  >cr 


mn 


<r 


mn 


(2-r»u) 


This  expression 


can  be  evaluated  further  us inn  the  Debye  approximation 


03  =  VK 

where  v  is  an  average  of  the  transverse  and  lonritudinal  velocities. 
Por  t’ne  behavior  at  T  =  0  (where  n=0)  then  from  eq.  (2-65)  eq . 

(2-6*1)  may  be  evaluated 


-ft*  2  r 

„  ?.  =  i  v  __d_  r  1 

1  iJ  nv  L 


7~TP~ 

D  mn 


[oor’(KDrnn,-1]] 


For  hiqh  temperatures  n  -*kpT/fiw  so  that  eq.  (2-6*0 


k  Tk 

,2-1  2V_  __B _ D 

^  2n2  mv2 


1- 


S1(k  r  )  q 
D  mn 

<  r 
D  mn 


(2-66) 

becomes 

(2-67) 


where  Si  is  the  usual  sine  integral. 

The  second  term  in  the  brackets  in  both  eq.  (2-67)  and 

eq.  (2-67)  comes  from  the  cross  term  in  eq.  (2-5*0  and  both  no  to 

21  22 

zero  as  r-***.  These  results  aqree  with  Schmidt  5  in  the  two 
temperature  limits.  Schmidt  extends  the  calculations  for  the 
intermediate  temperatures  of  interest  for  EXAFS. 


MKMMSBMKJ 
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He  proposes  that  the  EXAPS  temperature 
shell  be  of  the  form  exp  (-1  )  where 


correction  for  the 


I 

.1 


=  i^L-k2 

kBm6D 


<P  (T/e  ) 
J  D 


(2-68) 


and  kB  is  the  Boltzman  constant,  m  is  the  mass  of  the  atoms  in  the 
absorber,  0D  is  the  Debye  temperature  of  the  absorber,  k  is  the 
wave  vector  of  the  photoelectron,  and 

VT/V  ‘  *2J  -|+  <2-63) 

wnere  4-^,  and  <j>3j  are  defined  in  reference  22  and  depend  only 

on  the  crystal  structure  of  the  absorber.  All  data  presented  herein 
were  taken  at  liouid  nitrogen  temoerature  (T  =  77.4°K).  Por  these 
materials  T/eD  x.2  while  Shmidt  sup;nested  that  en.  (2-68)  is  valid 
for  T/od  '''.33.  However,  there  is  little  difference  in  the 
calculated  EXAPS  between  the  correction  from  eq.  (2-68)  or  the 
result  derived  in  ref.  22  which  is  valid  for  T  «  e  so  that  eq. 

(2-68)  was  used  as  it  was  somewhat  easier  to  evaluate.  It  should 
be  noted  that  the  temperature  factor  from  en.  (2-68)  cast  in  a 

general  Debye  form  that  can  be  used  in  eq.  (2-34)  by  makin.r  the 
identification 


2_  6 


o .  =  - 


,  41 .  (T/e  ) 

kBraeD  3  D 


(2-70) 
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so  that  the  differences  in  coupling  coefficients  are  related  to 

the  differences  in  4>  for  each  shell.  In  materials  such  as 

amorphous  semiconductors  where  there  is  no  method  to  calculate  the 

thermal  smearing  for  each  shell,  it  would  be  possible  in  theory  to 

measure  EXAFS  as  a  function  of  temperature  and  fit  the  variations 

to  a  Debye  function.  Exnerimentall.v  this  has  been  done  for  some 

23,24,25,2 6 

crystalline  materials  ,  but  it  is  usually  difficult  to 

separate  out  the  shell  to  shell  differences.  Thus,  the  attempts 
at  experimentally  determining  the  temperature  variation  have  met 
with  little  success.  In  the  next  chapter  it  is  shown  how  structural 
information  can  be  obtained  directly  from  EXAFS,  and  it  is  oUg— 
gested  how  an  experimental  measurement  of  the  temperature  variations 
will  be  possible.  Finally,  in  the  amorphous  materials  there  is 
also  a  contribution  to  the  average  displacement  of  an  atom  about 
its  equilibrium  position  because  of  structural  disorder.  Thus  a 
has  both  static  and  dynamic  contributions  which  must  be  separated 
in  any  structural  analysis.  This  will  also  be  discussed  in 
Chapter  3- 

As  mentioned  the  phase  shift,  ri(k),was  calculated  using  the 

WKB  method  from  eq.  (2-24).  The  self-consistent  Hartree  atomic 

27 

potentials  of  Herman  and  Skilman  were  used  for  all  calculations. 
These  atomic  potentials  will  be  shielded  in  the  solid  and  will 
not  have  as  large  a  tail  as  that  given  by  Herman  and  Skilman,  thus 
providing  an  error  in  n  .  An  additional  error  is  introduced  by  the 
fact  that  the  potential  of  the  ionized  atom  differs  from  that  given 
by  Herman  and  Skilman  because  of  the  missing  K  shell  electron. 
Further  error  is  introduced  by  the  WKB  approximation.  Eor  these 
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reasons  the  calculations  of  n  are  not  expected  to  be  accurate 

enough  to  give  detailed  agreement  with  the  experiment,  although 

one  expects  semi-quantitative  agreement.  We  use  the  WKB  method 

as  a  tentative  example  of  phase  shift  determination.  The  final 

values  of  n  were  determined  from  the  experimental  data  and  then 

used  to  calculate  the  fine  structure.  This  was  done  as  follows: 

i he  primary  peak  positions  were  assumed  to  come  from  the  first 

coordination  shell.  The  energy  positions  of  these  peaks  were 

usea  to  calculate  the  phase  shifts  for  these  energies  using  eg. 

(2-34).  This  was  done  for  various  choices  of  T  in  en.  (2-1 6) 

KK  * 

and  the  results  compared  to  the  Phase  shift  curve  determined  from 
the  WKB  approximation.  The  WKB  curve  and  the  experimentally 
determined  curve  for  germanium  are  shown  in  Pig.  3.  This 
technique  is  similar  to  methods  used  by  Kozlenkov9  and  by  fiott17, 
and  it  underlines  the  close  relationship  between  the  atomic  potential 
the  inner  crystalline  potential,  and  the  phase  shift.  If  the 
nroper  potentials  aren't  known,  then  k  and  n  must  be  determined 
self-consistent l.v  by  comparing  the  result  of  eq.  (2-34)  to  the 
experimental  EXAFS.  1/e  also  emphasize  that  more  suitable  trial 
potentials  than  from  Herman  and  Skilman  may  be  found  for  some 
materials ;  however,  the  Herman  and  Skilman  potentials  represent  a 
convenient  starting  point  for  all  atomic  species  studies  by  EXAFS, 
especially  since  the  phase  shift  is  most  sensitive  to  the  potential 
of  the  type  of  absorbing  atom  and  not  on  its  crystal  structure, 
bonding  state,  or  kinds  of  neighbors  (e.g.,  the  phase  shift  o/ne 
in  crystalline  Oe,  OeOg  and  OeSe  are  found  to  be  the  same  to  within 
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Figure  3 

Comparison  of  the  WKB  phase  shift  using  the  Herman-Ski liman 
potential  and  the  fitted  linear  phase  shift  for  Ge. 


experimental  error) . 

Theoretical  EXAFS  was  calculated  using  eq.  (2-34)  for  Cu 
(f.c.c.  lattice),  a-Fe  (b.c.c.)  and  fle  (dia.  cubic)  and  compared 
to  experiment.  The  data  used  for  each  element  are  listed  in 
'able  II  where  W Q-\J  from  eq.  (2-35),  nOO,  e^,  and  the  constant 
factor  In  the  temperature  correction  (eq.  2-55)  are  listed  along 
with  the  values  of  N,  r,  and  a  for  thirteen  coordination  shells. 

The  phase  shifts  are  linear  of  the  form 

n  =  -ak  +  b  (2-71) 

It  has  been  found  that  summing  to  thirteen  coordination  shells  is 
sufficient  to  Include  all  shells  which  give  a  sigrtificant  contri¬ 
bution  to  the  fine  structure.  The  summation  of  eq.  (2-34)  has 
also  been  performed  to  Include  20  shells  and  no  significant  change 
in  EXAFS  was  observed. 

In  Fig.  4a  and  b  comparison  of  the  theoretical  and  experi¬ 
mental  EXAFS  of  Cu  and  Fe  Is  shown.  The  theoretical  curves  are 
calculated  using  eq.  (2-34)  and  the  data  of  Table  II.  The  experi¬ 
mental  curves  were  from  well-annealed  foils  of  2.5p  and  6P  thickness 
for  Cu  and  Fe  respectively.  The  data  were  all  taken  at  77°K. 

These  curves  are  the  oscillatory  part  of  the  absorption  coefficient, 
corresponding  to  the  definition  of  x  in  eq.  (2-34).  The  sloping 
background  of  the  absorption  coefficient  has  been  subtracted 
graphically.  A  comparison  of  the  curves  shows  very  good  agree¬ 
ment  In  the  position,  shape  and  relative  heights  of  the  EXAFS 
features  over  the  entire  experimental  range.  Belov;  100  eV  there 
Is  a  disagreement  in  the  heights  of  the  features  because  the 
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Table  II 


Input  Data  to  Calculate  EXAPS  in  Cu,  Pe,  and  Oe 


Cu 

fee 


3.67ft 


k  = 

263(E+10)a 

n  = 

-.245k  +  7.95 

y  = 

.3 

e  = 

320°K 

D  b 

12fl2/kBn9D  = 

.02643 

b 

Shell,  ,j  = 

N  r  /a  <j> . 

J  3  J 

1 

12  0.707  .2333 

2 

6  1.000  .2800 

3 

24  1.225  .2809 

n 

12  .1414  .2797 

5 

24  1.581  .2838 

6 

8  1.732  .2904 

7 

48  1.871  .2960 

8 

6  2.000  .2992 

9 

36  2.121  .3003 

10 

24  2.236  .3002 

11 

24  2.345  .2998 

12 

24  2.449  .2998 

13 

72  2.550  .3005 

a)  The  zero  of  energy  is  taken 
inition  of  k  is  taken  from 


Fe 

Oe 

bee 

dia.  cubic 

2.85ft 

5.65ft 

.263(E+15)a 

.263(E-l)a 

-.194k  + 

7.96 

-.267k  +  7-60 

.3 

.3 

4l0°K 

370°K 

.02515 

.02149 

b 

b 

M  r  /a 

* 

N  r  /a  <J>. 

3  3 

3 

3  3  J 

8  0.866 

.2098 

4  0.433  .1694 

6  1.000 

.2355 

12  0.707  .2624 

12  1.414 

.2598 

12  0.829  .2723 

24  I.658 

.2570 

6  1.000  .2708 

8  1.732 

.2566 

12  1.096  .2698 

6  2.000 

.2601 

24  1.225  .2724 

24  2.208 

.2664 

16  1.299  .2758 

24  2.236 

.2672 

12  1.414  .2819 

24  2.249 

.2676 

24  1.479  .2850 

32  2.598 

.2740 

24  1.581  .2882 

12  2.828 

.2740 

12  1.639  .2890 

48  2.958 

.2736 

8  1.732  .2893 

30  3.000 

.2735 

24  1.785  .2890 

to  be  the 

Pe  rmi 

energy.  This  def- 

eq.  C 2 — 1 6 ) 

and  (2 

-35) 

b)  Prom  ref.  21  and  22  and  eq.  (2-68)  and  (2-69) 
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Figure  4 


Comparison  of  theoretical  and  experimental  EXAFS  for 
copper  and  iron.  The  theoretical  spectrum  is  obtained 
from  eq.  (2-34)  and  the  data  of  Table  II. 
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PHOTOELECTRON  ENERGY,  eV 


FI pure  5 

Comparison  of  theoretical  and  experimental  EXA^S  for  permanium. 
Curve  (a)  is  the  experimental  data  with  the  usual 
photoelectric  edpe  removed,  curve  (b)  is  the  smoothed 
experimental  curve,  and  curve  (c)  is  the  theoretical  curve 
calculated  from  eq.  (2-3*0  and  the  data  of  Table  II. 
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Figure  6 

(a)  Theoretical  EXAFS  for  iron  without  the  temperature  and 
ran^e  terms,  (b)  temperature  correction  has  been  apDlied 
to  curve  (a),  (c)  the  exponential  attenuation  has  been  applied 
yielding  the  theoretical  curve  in  Fig.  4b.  All  curves  are 
plotted  to  the  same  (arbitrary)  scale. 
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scattering  factor  used  in  eq.  (2-34)  was  not  adequate.  There  is 
also  some  disagreement  in  shape,  at  the  lower  energies  around  50 
eV,  but  as  mentioned  previously  the  approximations  made  were  not 
expected  to  be  valid  in  this  range. 

In  Fig.  5  a  comparison  of  the  theoretical  and  experimental 
EXAFS  for  germanium  is  shown.  Curve  (a)  is  the  oscillatory  part 
of  the  experimental  absorption  coefficient  including  statistical 
noise.  Curve  (b)  is  a  smoothed  version  of  Curve  (a)  where  the 
statistical  noise  was  filtered  out  using  a  computer  method  which 
will  be  described  in  Appendix  A,  and  Curve  (c)  is  the  theoretical 
curve  calculated  from  eq.  (2-34)  using  the  data  of  Table  II.  The 
agreement  between  the  theoretical  and  experimental  spectra  is 
excellent  over  the  entire  range.  Even  below  100  eV  there  is 
general  agreement  with  the  position  of  the  maxima,  although  the 
heights  of  the  peaks  In  the  theoretical  spectrum  are  too  large, 
which  as  mentioned  above,  may  be  due  to  the  inadequate  scattering 
factor  used  at  the  lower  energies,  the  neglect  of  band  effects 
on  the  density  of  states,  or  that  the  experimental  resolution 
of  4-6  eV  was  not  included  in  the  calculated  spectrum.  This  . 
would  have  the  effect  of  smoothing  and  reducing  the  features  and 
would  be  most  effective  at  the  lower  energies  where  the  widths 
of  the  features  are  only  a  few  eV  wide. 

In  Fig.  6  the  importance  of  the  temperature  and  range  terms 
in  calculating  EXAFS  which  agree  with  the  experimental  data  is 
shown.  Curve  (a)  is  the  curve  for  iron  from  Fig.  4b  without  the 
temperature  or  the  range  term.  Curve  (b)  is  the  same  as  Curve  (a) 
except  the  temperature  factor  from  eq.  (2-68)  is  applied.  In 


Curve  (c)  the  range  term  is  also  applied,  yielding  the  theoretical 

curve  of  Pig.  4b.  The  importance  of  both  corrections  in  lowering 

the  amplitudes  of  the  peaks  especially  at  higher  energies  and 

in  smoothing  the  spectrum  by  reducing  the  contributions  of  higher 

coordination  shells  can  be  easily  seen.  It  is  obvious  why  earlier 
3  9 

attempts  *  to  make  quantitative  comparisons  of  the  theoretical 
and  experimental  data  without  applying  these  corrections  met  only 
limited  succe-s. 

The  computer  nrogram  used  to  calculate  theoretical  EXAPS  from 
eo.  (2-34)  is  listed  in  Aopendix  C. 


Chapter  3 

A  New  .Method  for  Structure  Analysis  Using; 
The  Fourier  Transform  of  EXAFS 


In  Chapter  2  a  simple,  flexible  theory  of  EXAFS  was  presented 
which  allows  calculation  of  fine  structure  spectra  which  are  In 
good  agreement  with  observed  spectra.  Perhaps  of  more  practical 
interest  is  the  inverse  problem  of  determining  the  atomic  arrange¬ 
ments  in  the  material  from  its  EXAFS;  this  would  constitute  a 
new  and  valuable  research  tool.  The  results  of  the  point  scattering 
theory  (eq.  (2-3*0)  showed  that  EXAFS  was  given  by  a  sum  of 
damped  sine  waves  each  originating  from  a  coordination  shell. 

This  suggested  that  the  data  could  be  inverted  using  Fourier 
analysis . 

If  eqX  2-3*0  is  rewritten 


-x(k) 


-1  ,  -1  V  j 

k  f  (k)  = 


-2oj2kZ 


sin  ^  2kr  +  2n  (k)j 
J 


(3-D 


the  Fourier  transform  of  this  equation  is 


f x(k)k”1f(k)"1  sin(2kr+2n)dk  = 

o  (r-r.)2 


I 


U,e-YrJ 

8r  ^  o 

5  J 


2 


e 


+  A  (r)  =  <t>(r) 


(3-2) 
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Here  A(r)  is  a  function  which  is  sharply  peaked  about  -r  and  is 
very  small  in  the  space  r  >0  of  physical  interest  and  will  be 
neglected  from  now  on.  The  right  hand  side  of  eq.  (3-2),  <Kr), 

is  a  radial  structure  function.  It  shows  that  the  transform  of 
the  experimental  data,  given  by  the  left  hand  side  of  eq.  (3-2), 
is  a  series  of  normalized  Gaussian  curves  whose  amplitudes  are 


8r  2o 
J  J 


(3-3) 


and  whose  widths  W 

j 


value  to  A  /e  are 


where  the  amolltude  is  down  from  its  peak 


W  =2/2  a 

J  j 


(3-4) 


Note  that  in  crystalline  materials  the  atoms  in  a  coordination 
shell  will  be  exactly  at  r j  and  the  ,ys  determined  from  the  data 
using  eq.  (3-4)  are  entirely  thermal  in  origin.  For  simple 
materials  such  as  Fe,  Cu  and  Ge,  the  measured  o  '  s  can  be  compared 
to  the  <ys  calculated  from  Shmidt's  theory21-2^.  In  more  oonlpl(,x 
materials  the  widths  of  the  Gaussians  can  be  measured  for  each 
shell  as  a  function  of  temperature,  and  fitted  to  Shmidt's  theory 
or  a  similar  function.  This  could  be  used  to  determine  a  Debye 
temperature  for  the  material. 

For  structure  analysis  in  crystalline  materials  the  0  's 
can  be  measured  or  calculated  and,  since  the  N  »s  ancj  r  »s 

known  from  the  crystal  structure,  the  range  parameter  y^can  be 
determined  from  the  measured  amplitudes  A  using  eq.  (3-3).  This 
can  be  seen  by  rearranging  eq.  (3-3),  ^ 
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-In  A  /(N  /8r  2o  )  =  yr  (3-5) 

J  j  J  J  J 

Thus,  plotting  the  left  side  of  eq.  (3-5)  vs.  r  yields  a  straight 
line  whose  slope  is  y  .  This  will  be  done  in  Chapter  5  for  Ge. 

It  must  be  emphasized  that  this  method  is  the  only  reliable  way 
to  determine  the  value  of  y.  Other  calculations  of  y  are  only 
approximate . 

Another  parameter  that  can  be  obtained  from  the  radial  struc- 

1 

ture  function  of  crystalline  materials  is  the  phase  shift.  In 
deriving  eq.  (3-2),  it  was  assumed  that  the  nhase  shift  was  known 
and  included  in  the  argument  of  the  analyzing  waves.  Often, 
however,  the  phase  shift  is  not  known  sufficiently  well  to  do  this; 
thus,  the  argument  of  the  analyzing  wave  used  is  2kr  instead  of 
(2kr  +  2n).  The  result  of  the  theory  given  by  the  right  hand  side 
of  en.  (3-1)  contains  the  phase  shift  which  is  assumed  to  have  the 
linear  form  given  by  eq.  (2-71)  so  that  when  the  Fourier  integra- 

■] 

l 

tlon  is  done  the  right  hand  side  of  eq.(3-2)  remains  the  same 

2 

except  that  the  argument  of  the  exponential  is  now^r-(r  -a)^  , 

J 

where  a  is  the  k  dependent  part  of  the  phase  shift.  This  has  the 
effect  of  shifting  all  of  the  peaks  toward.  the  origin  by  an 
amount  a,  where  a  is  usually  less  than  .4  so  that  a  <<  r  in  all 
cases.  From  the  known  distances  of  the  coordination  shells  in  a 

f 

given  crystal  structure,  the  amount  of  the  shift  can  be  deter¬ 
mined  for  that  particular  element.  Once  the  k  dependence  of  the 
phase  shift  is  known,  the  constant  part  of  the  phase  shift  can 
be  determined  from  the  experimental  data.  Depending  on  how  the 
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transform  is  taken  it  may  not  be  necessary  to  know  b  to  invert 
the  data  however,  it  is  convenient  to  have  to  generate  the  linear 
phase  shift  so  that  EXAPS  can  be  calculated  using  eq.  (2-34)  or 
the  phase  shifts  for  the  same  kinds  of  atoms  in  different  material 
can  be  compared.  Finding  b  is  done  by  locating  the  extrema  of  the 
peaks  due  to  the  dominant  shell  (usually  the  first  shell)  and 
using  the  equation  for  the  extrema  of  the  sine  wave 


2k  (r  -a)  +  2b  =  (2n  +  1)  -  (3-6) 

n  j  2 

where  a  and  be  are  the  phase  shift  parameters  from  eq.  (2-71), 

t  hi 

kn  is  the  k  Value  of  the  n  extremum,  and  n  =  0,2,4,...  for 
maxima  and  1,3,5,...  for  minima.  Since  the  absolute  value  of  n 
isn't  known,  eq.  (3-6)  can  be  solved  for  b  using  successive  pairs 
of  extrema  and  a  "best"  b  is  determined  from  the  result. 

In  noncrystalline  materials  the  problem  of  determining 
information  is  much  more  difficult  since  there  is  no  reliable  way 
of  measuring  or  calculating  parameters  such  as  n,  y,  and  0  in 
these  materials.  Furthermore  the  peak  width  a  contains  both 
static  (disorder)  and  dynamic  (thermal)  broadening.  If  Gaussian 
shapes  are  assumed  for  both  contributions  to  the  width,  then 


TOT2 

o 

J 


0  +  0 


J  j 


(3-7) 


where  T  and  D  stand  for  thermal  and  disorder  respectively.  Hence, 
for  a  given  shell  o^,  and  the  kinds  of  atoms  in  a  shell  are  all 

the  important  structural  information.  One  method  to  obtain 
structural  information  from  non-crystalline  materials  is  to 


compare  the  non-crystalline  radial  structure  function  to  the 
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function  from  its  crystalline  polymorph.  If  the  short  range  order 
is  similar  in  the  two  materials,  then  the  assumption  is  made  that 
the  range  parameter  y  is  the  same  in  both  the  crystalline  and 
non-crystalline  material.  This  should  be  true  if  the  near  neighbor 
environment  is  similar  in  both  materials  and  1/Y ,  the  mean  free 
path,  is  on  the  order  of  the  inter-atomic  spacing.  It  is  also 
assumed  that  thermal  smearing  for  similar  shells  is  the  same 
because  of  the  similarity  of  the  bonding  involved.  For  any  new 
peaks  appearing  in  the  amorphous  structure  function  this  assumption 
would  be  questionable,  but  this  would  be  most  serious  only  for 
the  first  one  or  two  shells  since  beyond  that  the  shell  to  shell 
variation  in  a  is  only  a  few  percent.  Finally,  It  is  assumed  that 
the  phase  shift  for  the  two  polymorphs  Is  the  same  since  it  has 
been  found  to  depend  primarily  on  the  potential  of  the  absorbing 
atom  and  not  on  the  surrounding  structure. 

These  assumptions  make  it  possible  to  obtain  structural 
information  from  the  measured  height  and  width  using  eqs.  (3-3) 
and  (3-4)  with  the  crystalline  values  of  y  and  o“,  and  the  fol- 
lowing  argument  shows  that  absolute  numbers  of  atoms  In  each  shell 
may  be  determined  from  the  data.  Even  though  phe  theory  does 
not  predict  an  absolute  value  of  the  EXAFS  cross  section,  it  is 
expected  that  it  is  In  error  only  by  an  overall  multiplicative 
constant  since  all  of  the  functional  dependence  was  included  in 
the  theory.  One  method  to  obtain  numbers  of  atoms  directly  from 
the  data  is  to  normalize  the  observed  jump  in  both  the  crystalline 
and  amorphous  materials  of  the  absorption  coefficient  at  the  K 
edge  to  the  same  (arbitrary)  value.  This  In  effect  normalizes  the 


cross  section  for  absorption  in  each  material,  so  that  the  overall 
constant  for  both  the  crystalline  and  non-crystalline  materials 
should  be  the  same.  If  the  cross  sections  have  been  normalized, 
then  the  results  of  the  theory  state  that  the  relative  amplitude 
of  the  observed  fine  structure  from  a  given  shell  will  be 
proportional  to  the  number  of  atoms  in  that  shell  in  the  two 
materials.  To  get  a  quantitative  comparison  of  number  of  atoms  in 
a  shell,  the  transforms  of  the  two  materials  should  be  compared 
using  eq.  (3-3).  The  overall  constant  missing  from  eq.  (3-3) 
should  be  the  same  and  so  the  absolute  values  of  the  A  '  s  obtained 
from  the  data  can  be  used  to  obtain  the  number  of  atoms  in  an 

amorphous  shell  if  the  number  of  atoms  in  the  comparable  crystalline 
shell  is  known,  using 


N  (A)  =  N  (C) 

J  J 


A  (A)o  (C)  r  (C)2  -yr, (A) 
J  _ J  e  J 

(C)o  (A)  r  (A)2  e“YrJ(C) 


(3-8) 


In  many  amorphous  materials  the  first  shell  is  very  similar  to 

the  first  crystalline  shell  because  the  near  neighbor  bonding  is 

the  same,  so  that  r  (A)  =  r  (C)  and  a  (A)  =  0  (c).  Therefore  a 

J  J  J  j 

simplified  form  of  eq.  (3-8)  may  be  used  to  find  N  (A) 

J 


51 


Once  N  Is  known  for  any  shell  In  the  amorphous  material,  all  of  the 

j 

other  N  's  can  be  determined  by  comparing  the  amplitude  of  their 
shells  to  the  shell  where  N  is  known  and  applying  eq.  (3-3).  Thus, 
even  though  absolute  absorption  data  isn't  taken,  the  EXAPS 
technique  provides  a  straightforward  method  of  determining  the 
number  and  distribution  of  atoms  about  a  given  distance  in  most 
amorphous  materials. 

The  preceding  discussion  of  the  kinds  of  structural  information 

obtainable  from  EXAPS  data  was  based  in  part  on  the  fact  that  the 

integration  in  eq.  (3-2)  was  over  all  k  space.  In  fact,  the 

EXAFS  data  is  taken  only  over  a  small  domain  in  k  space  which 

typically  is  from  k  ^2,  because  the  data  near  the  edge  often 

contains  features  from  effects  other  than  EXAFS,  to  k  ^20,  which 

is  the  practical  experimental  limit,  since  the  time  involved  to  go 

farther  would  be  prohibitive.  When  a  Fourier  transform  is  made 

over  a  finite  range  a  termination  error  is  introduced  into  the 

radial  structure  function.  The  corresponding  problem  in  radial 

28 

distribution  analysis  has  been  discussed  by  Warren  who  shows  that 
this  error  reduces  the  heights  of  the  peaks,  broadens  them  some¬ 
what,  and  introduces  a  set  of  ripples  on  either  side  of  the  peak. 
Because  the  precise  effect  of  the  termination  error  in  radial 

distribution  analysis  is  given  by  a  complicated  convolution 
28 

integral  ,  the  data  is  usually  corrected  for  termination  error 

by  either  applying  an  exponential  convergence  factor  which  reduces 
28  “ 

the  rinples  ,  or  by  a  complicated  iterative  procedure  designed  to 

29 

experimentally  subtract  the  termination  ripple  from  the  data  , 
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or  by  using  analytic  continuation  of  the  data3  to  extend  the 
range  of  the  data  and  thereby  reduce  the  termination  effects. 

The  effect  of  termination  error  on  the  transform  of  EXAFS 
data  is  shown  in  Fig.  7,  where  the  theoretical  He  EXAFS  spectra 
of  Fig.  5c  has  been  transformed  using  eq.  (3-2),  except  that  the 
upper  limit  of  integration  is  =  15.  As  is  evident  from  the 
figure  there  is  some  ripple  introduced  about  each  peak  as  a  result 

of  the  transform. 

One  difference  in  the  EXAFS  termination  error  from  that  of 
radial  distribution  analysis  is  that  the  effect  of  termination  can 
be  evaluated  by  direct  integration.  The  integral  of  interest  is 

X  =  1  fa  e  ^  cos  2kAr  dk  (3-10) 

v/2n  j  J 

o 

where  a  =JTo  ,  Ar  =  r-r  ,  and  k  is  the  maximum  value  of  k  as 
j  1  j  j  m 

determined  by  the  range  of  data  used  in  evaluating  the  integral 
from  the  left  hand  side  of  eq.  (3-2).  The  int&gral  (3-10)  is  the 
term  giving  the  series  of  (laussians  in  the  region  r  >0.  The  other 
term  which  would  have  given  peaks  about  r  —r^  may  still  be  ignored 
since,  even  with  termination  effects,  it  will  have  no  appreciable 
effect  for  r  >0.  The  integrand  from  (3-10)  may  be  rewritten  using 
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Figure  7 


Termination  effects  on  the  radial  structure  function  of  the 
theoretical  Ge  EXAFS  of  Fip;.  5c.  (k  *  15). 
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where  the  same  k  independent  term  used  to  complete  the  square  has 

been  factored  from  both  integrals.  These  integrals  may  be  evaluated 

in  terms  of  error  functions  after  changing  the  variables  to  = 

a  k  +  iAr  /a  and  x  =  a  k  -  iAr  /a 
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where  the  error  functions  from  the  bottom  terms  have  cancelled 

since  erf(-x)  =  -erf(x).  This  expression  may  be  reduced  even 

further  >b,y  using  the  expansion  of  an  error  function  of  complex 

30 

argument  given  in  Abramowitz  and  Stegun  f  So  that 
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Although  this  expression  is  complicated,  several  factors  are 
important.  First,  unlike  radial  distribution  analysis,  where  the 
effect  of  termination  is  described  by  a  convolution  integral, 
here,  the  termination  error  merely  multiplies  the  Gaussian  curve. 
Secondly,  the  peak  amplitude  (at  Ar  -  0)  is  reduced  by  a  factor 

J 

of  erf (a  k  ),  where  typically  a.k  *1.  Away  from  the  peak  the 
J  m  J  m 

termination  introduces  a  ripple  of  periodicity,  n/k^. 

A  comparison  of  the  termination  error  that  is  calculated 
using  eq.  (3-16)  versus  the  error  obtained  from  direct  integration 
is  shown  in  Fig.  8.  The  parameters  for  the  original  Gaussian  are 
shown  in  the  figure  and  correspond  to  the  first  shell  of  Ge.  The 
expression  containing  termination  ripple  calculated  from  eq.  (3-15) 
is  shown  in  curve  b  and  the  first  shell  curve  obtained  from  direct 
integration  as  in  Fig.  7  is  shown  in  curve  c.  The  reduction  in 
amplitude  and  the  broadening  of  the  peak  expected  from  eq.  (3-15) 
are  shown  in  the  figure  where  a  comparison  of  curves  b  and  c  show 
good  agreement  between  the  calculated  and  integrated  result.  The 


Figure  0 

(a)  First  shell  (laussian  from  the  theoretical  he  rsf, 

(b)  calculated  termination  effect  using  eq.  (2-15)  for  k  = 

(c)  termination  effects  on  (a)  from  direct  integration  with  k 
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differences  between  b  and  c  above  .  25A  are  the  result  of  an 
approximation  used  in  evaluating  eq.  (3-15)  in  this  region.  Unfor¬ 
tunately,  the  presence  of  statistical  noise  in  the  data  and  inter¬ 
polation  between  data  points  during  the  integration  routine  has 
made  it  impossible  to  use  eq.  (3-15)  to  correct  the  actual  experi¬ 
mental  data  for  termination  error.  The  results  of  this  equation 
can  be  used  to  obtain  the  height  and  width  of  the  Gaussians  without 
termination  for  the  major  peaks,  but  it  cannot  be  used  to  identify 
or  correct  small  peaks  whose  amplitudes  are  on  the  order  of  the 
noise.  These  peaks  are  important  in  the  study  of  non-crystalline 
materials  and  future  work  should  develop  techniques  which  will 
remove  or  reduce  the  termination  effects  and  allow  the  resolution 
of  small  signal  peaks  from  the  data. 

One  final  point  should  be  made  about  comparing  the  data 
obtained  from  the  EXAPS  technique  versus  the  data  obtained  from 
radial  distribution  analysis.  The  resolution  possible  in  any 
transform  atom  of  scattering  or  EXAPS  data  depends  on  the  number 
of  oscillations  the  various  shells  contribute  over  the  range  of 
the  data.  If  EXAPS  results  are  compared  to  any  other  method  of 
determining  radial  distribution  functions  (x-ray,  electron,  or 
neutron  scattering),  and  the  data  are  taken  over  the  same  range  of 
k  space,  the  scattering  intensity  data  varies  as  sin  kr,  and  the 
EXAPS  data  varies  as  sin  2kr.  Thus,  there  are  twice  as  many 
oscillations  in  the  EXAPS  data,  hence  twice  the  resolution  in  the 
transformed  data  will  be  possible. 


Chapter  4 

Experimental  Methods 


In  this  chapter  the  experimental  methods  for  obtaining 
EXAPS  data  are  described.  Although  many  techniques  have  been 
used  throughout  the  years  to  record  fine  structure  spectra  the 
results  have  generally  been  of  limited  value  because  of  the 
poor  resolution  and  counting  statistics,  limited  energy  range 
above  the  absorption  edge,  and/or  very  slow  data  collection. 

The  system  described  here  permitted  automated  collection  of 
medium  resolution  data  with  good  counting  statistics  over  energy 
intervals  up  to  2  keV  above  the  absorption  edge.  All  of  this 
was  necessary  to  obtain  useful  structural  information  from  the 
Fourier  transform  of  the  data. 

The  basic  x-ray  system  was  a  Siemen's  Crystalloflex  4 
power  supply  and  horizontal  single-crystal  diffractometer.  The 
x-ray  continuum  radiation  from  medium  to  heavy  target  tubes, 
such  as  W,Ag,  and  Mo,  were  operated  at  their  highest  power 
levels  (usually  Ho  kV  and  30  mA)  for  maximum  intensity.  To 
achieve  adequate  resolution,  precise  alignment  of  the  spectro¬ 
meter  was  obtained  using  narrow  slits  (usually  50  u).  The 
analyzing  crystal  used  was  a  carefully  selected  single  crystal 
of  kiF  (2d20Q  =  0265A)  mounted  in  a  holder  which  allowed 

accurate  placement  and  alignment  of  the  crystal  in  the  beam. 

The  x-rays  were  detected  by  a  Nal  (Tl)  scintillation  counter 
with  a  pulse  height  analyzer  to  eliminate  higher  order  reflect¬ 
ions,  and  recorded  on  a  scaler/timer  operated  in  a  fixed  count 
mode  to  assure  good  counting  statistics  for  each  data  point. 
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Usually  the  time  to  accumulate  105  counts  ('••3*  statistical 
counting  error)  was  measured.  The  numerical  resolving  power 
of  the  spectrometer,  which  corresponds  to  an  energy  resolution 

of  4-6  eV  was  A/AX  ^  2000  -  3000. 

The  sample  was  mounted  in  a  cryostat  and  suspended  In  the 
beam  between  the  entrance  slit  and  analyzing  crystal  on  a  holder 
which  can  automatically  be  moved  in  and  out  of  the  x-ray  beam.. 

The  technique  for  doing  this  is  illustrated  in  fig.  9.  The 
sample  holder  was  spring  loaded  to  return  to  one  position  and 
actuated  by  means  of  a  cam  in  contact  with  the  center  rod.  The 
total  excursion  was  controlled  by  mechanical  stops  with  the 
sample  holder  isolated  from  the  cryostat  by  a  flexible  bellows. 
The  cam  was  activated  by  a  connecting  rod  connected  about  1" 
off  center  on  the  shaft  of  a  Slo-Syn  motor  in  a  bell-crank 
arrangement.  The  rod  was  alternately  pulled  and  released  by 
rotating  the  shaft  180°  each  time  the  sample  was  to  be  moved. 

As  the  Slo-Syn  shaft  was  rotated  in  the  same  direction  each  time, 
a  separate  system  of  cams  ar.d  microswitch  actuated  a  stepping 
motor  driving  the  x-ray  goniometer.  The  goniometer  could  be 
stepped  in  Increments  of  .0025  -  .025° 20  depending  upon  the 
requirements  of  the  experiment.  Normal  step  sizes  were  .005° 
or  .010°.  A  schematic  diagram  of  the  apparatus  is  shown  in 

fig.  10. 

During  a  typical  experimental  run  the  goniometer  was 
positioned  at  a  fixed  starting  angle  and  the  absorber  positioned 
out  of  the  beam.  When  10^  counts  had  been  accumulated  a  signal 
was  sent  from  the  scalar  to  a  Hamner  Scanner-Coupler  which  held 
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Figure  9 


Diagram  of  the  sample  channinp;  mechanism. 
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Figure  10 


Schematic  of  the  experimental  apparatus, 
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the  counting  sequence,  signaled  a  Humphrey  Stepper  Controller 
to  rotate  the  Slo-Syn  motor  repositioning  the  sample,  and 
relayed  the  digital  information  to  an  IBM  keypunch  which  punches 
the  accumulated  time.  This  process  was  then  repeated  with  the 
absorber  in  the  beam.  Each  time  the  sample  changer  returned  to 
its  original  position  the  goniometer  was  stepped  by  .005°. 

This  process  was  automatically  repeated  until  the  entire  region 
above  and  through  the  absorption  edge  was  recorded.  Typically, 
a  run  extended  from  100  eV  below  the  edge  to  1500-2000 
above  the  absorption  edge.  This  corresponded  to  about  a  4°20 
range  in  the  energy  regions  usually  studied. 

In  order  to  obtain  structural  information  using  the  method 
described  in  Chapter  3  the  thermal  smearing  of  the  data  must 
be  minimized,  which  was  done  by  taking  most  data  at  liquid 
nitrogen  temperature  (77°K).  Although  the  cryostat  was  capable 
of  operating  down  to  4°K,  the  small  reduction  in  the  thermal 
motion  over  that  at  77°K  was  not  worth  the  added  effort  and 
expense.  Another  requirement  for  adequate  data  was  that  he 
noise  level  be  reduced  so  that  signal  could  be  extracted  from 
the  data  as  far  from  the  edge  as  possible  to  obtain  greater 
resolution  and  lessen  Fourier  termination  effects.  The  major 
contribution  to  the  noise  was  due  to  the  statistical  counting 
fluctuations.  Although  10-'  counts  limits  the  statistical 
fluctuations  in  the  measured  absorption  cross-section  to  ^.3%, 
it  must  be  remembered  that  the  signal  of  interest  was  the  small 
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oscillation  on  top  of  the  usual  photoelectric  cross-section. 

105  counts  was  the  practical  limit  for  the  number  of  counts 
collected  for  each  point  because  of  long  term  drift  in  the 
experimental  apparatus.  Improved  statistics  were  achieved  by 
gathering  data  from  several  runs,  usually  10  or  more,  and 
averaging  them  together.  This  was  done  since  the  noise  was 
random  and  the  starting  points  of  the  data  were  well  aligned. 
This  process  was  time  consuming  since  each  run  took  about  1  day. 
However,  the  final  statistical  fluctuations  were  %.1%  which 
allowed  resolution  of  the  signal  over  a  few  hundred  eV  more 
range . 

The  data  cards  collected  from  each  run  were  analyzed  by 

a  series  of  computer  programs.  The  first  program  checked  for 

errors  that  may  have  occurred  during  the  experimental  run. 

Typically  the  errors  were  mechanical  and  caused  an  I  or  Iq 

point  to  be  missed.  To  keep  the  (I,  I  )  pairs  in  proper  order, 

values  for  the  missed  I  or  I  were  estimated  and  inserted. 

o 

This  was  not  a  serious  problem  occurring  no  more  than  2  or  3 
times  per  run. 

All  of  the  corrected  data  decks  from  a  given  sample  were 
then  combined  into  a  single  data  deck.  The  careful  alignment 
necessary  was  done  by  using  the  characteristic  tube  lines  and 
W  impurity  lines  present  in  the  data.  There  were  usually 
several  of  these  lines  in  the  data  which  allowed  alignment  of 
all  of  the  decks  to  ^.0025°  which  reflects  the  uncertainty  of 
the  Initial  positioning  of  the  spectrometer. 


The  combined  data  deck  was  then  analyzed  by  the  main 
program  which  computes  ln(Io/I)  and  the  energy  scales  from  the 
angular  data.  The  program  also  made  internal  corrections  of 
the  energy  scale  based  on  the  positions  of  the  impurity  lines 
whose  energies  were  known32  and  corrected  for  the  coincident 
losses  when  intense  impurity  lines  appeared  in  the  data.  This 
program  also  had  the  capability  to  Fourier  analyze  the  data 
using  the  fast  Fourier  transform  (fft)  described  in  Appendix  A. 
The  purpose  of  the  analysis  at  this  point  was  to  remove  the 
fine  structure  data  from  the  total  spectrum.  Filtering  out  the 
lowest  frequency  from  the  transformed  data  and  then  retrans¬ 
forming  the  data  accomplished  this  task.  This  method  of 
extracting  the  fine  structure  data  was  preferred  because  of  the 
simplicity  of  the  technique.  Since  the  shape  of  the  photo¬ 
electric  spectrum  was  not  accurately  known  near  the  edge,  a 
complicated  fitting  program  would  be  required  for  each  edge 
in  order  to  find  the  correct  background  spectrum  to  remove  from 
the  data.  The  fft  program  was  also  used  to  filter  out  much  of 
the  statistical  noise  in  the  data  to  produce  curves  as  in  fig. 
5b.  This  was  done  only  to  make  comparisons  easier;  when  the 
data  was  transformed  to  obtain  structural  information,  the 

statistical  noise  was  left  in. 

Proper  sample  preparation  was  also  very  important  to 
obtain  high  quality  EXAFS.  Samples  must  be  of  uniform  thick¬ 
ness  over  an  area  of  about  1cm  by  2.5cm.  The  metal  samples 
were  rolled  into  thin  foils  of  about  2-20y  thickness.  The 
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optimum  thickness  depended  on  the  absorption  coefficient  of 
the  material  but,  empirically,  it  was  observed  that  when  Iq/I 
3  just  above  the  absorption  edge,  the  best  resolution  was 
obtained.  For  alloys  and  compounds  the  situation  was  more 
complicated;  but,  in  general,  the  absorption  jump  was  made 
sufficiently  large  so  that  the  fine  structure,  whose  amplitude 
at  most  was  a  few  percent  of  the  absorption  jump,  could  be 
resolved.  The  rolled  foils  were  well  annealed  to  remove  any 
effects  of  work  hardening.  Other  materials  were  crushed  into 
a  fine  uniform  powder,  mixed  in  Duco  cement,  and  then  thin 
layers  were  spread  on  A1  foil.  (The  A1  foil  is  1  mil  thick 
and  used  as  a  backing  for  thermal  contact  with  the  sample 
holder,  but  had  a  low  absorption  coefficient  with  no  fine 
structure  in  the  energy  range  of  interest  (5-20  keV) . )  Care 
must  be  taken  to  make  uniform  samples  with  on  inhomogeneities 
or  pinholes.  Several  of  these  thin  layers  were  used  to  make 
uniform  samples  of  optimum  thickness.  Some  materials  which 
were  difficult  to  roll  or  grind  and  many  amorphous  samples 
were  evaporated  in  vacuum  onto  1  mil  A1  foil  or  mylar  to  a 
thickness  of  *1 u.  Multiple  layers  were  used  for  optimum 


thickness . 


Chapter  5 

Application  of  EXAFS  Technique  to  Structure  Determination 
of  Amorphous  Ge,  GeO^,  and  GeSe 

In  this  chapter  the  techniques  for  structure  determination 
using  EXAFS  will  be  applied  to  three  amorphous  materials — Ge, 

GeO^ ,  and  GeSe.  The  analysis  of  the  amorphous  samples  is  effected 
by  standardizing  the  EXAFS  parameters  with  their  respective 
crystalline  polymorphs,  as  described  earlier. 

The  techniques  to  reduce  the  effect  of  termination  errors 
have  not  yet  been  incorporated  in  the  data  analysis  programs; 
therefore,  the  information  from  the  peak  widths  is  unreliable. 
Considerable  structural  Information  can  be  derived  from  only  the 
positions  and  amplitudes  of  the  peaks. 

A.  Germanium 

Amorphous  Ge  was  the  first  system  studies  with  the  EXAFS 
technique.  It  has  been  studied  extensively  by  other  techniques, 
since  it  is  a  simple,  easily  prepared  amorphous  system.  The 
basic  structural  unit  is  the  tetrahedron  formed  from  Ge  making 
4  covalent  bonds.  In  crystalline  Ge  adjacent  tetrahedra  are  all 
linked  in  the  staggered  configuration  which  generates  the  diamond 
cubic  lattice.  In  amorphous  Ge  the  linking  of  the  tetrahedra  Is 
not  as  regular  leading  to  the  breakdown  of  long-rance  order. 

Two  primary  explanations  have  been  offered  to  describe  the  struc¬ 
ture  of  amorphous  Ge.  The  first  was  by  Grigorovici  and  Manaila33 * 3^ 
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who  pictured  the  structural  unit  as  a  pentagonal  dodecahedron 
made  of  12  five-fold  rings  of  linked  eclipsed  Ge  tetrahedra. 

Fig.  11  shows  the  staggered  and  eclipsed  linking  of  Ge  tetrahedra. 
Their  proposed  structure  consists  of  a  50-50  mixture  of  5  fold 
rings  and  distorted  crystalline  regions,  all  linked  together. 
Although  this  hypothetical  structure  is  in  general  agreement  with 
the  experimental  radial  distribution  function  (r.d.f.),  particu¬ 
larly  in  the  reduction  of  the  third  peak,  there  is  no  detailed 
agreement  and  the  density  appears  to  be  too  low.  Generally,  the 
structure  is  somewhat  artificial  and  the  requirement  that  only 
staggered  or  eclipsed  configurations  be  allowed  is  too  rigid  for 
a  random  structure. 

A  more  recent  model  has  been  proposed  by  Polk3^  based  on 
the  continuous  random  network  proposed  by  Zachariasen^^ .  Polk 
constructed  a  model  to  satisfy  local  bonding  requirements,  and 
minimize  strain  and  surface  density  of  unsatisified  bonds.  He 
concluded  that  his  Mo  atom  model  could  be  extended  indefinitely 
without  postulating  regions  of  special  bonding.,  The  r.d.f. 
calculated  from  his  model  agrees  very  well  with  the  r.d.f.  of 
amorphous  Si  by  Moss  and  Graczyk37  and  the  r.d.f.  of  amorphous 
Ge  by  Shevchik  which  show  the  first  and  second  peak  at  the 
crystalline  distances  and  a  small  third  neighbor  peak  at  a 
slightly  larger  distance  than  in  the  crystalline  case.  In  Polk’s 
model  the  tetrahedral  bond  angle  was  shown  to  deviate  1  10° 
about  the  ideal  angle  of  107°26'  and  the  rotations  of  the  tetra¬ 
hedra  about  the  common  bond  contained  intermediate  positions  in 
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Starve red  and  eclipsed  linking  of  He  tetrahedra. 
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addition  to  the  staggered  and  eclipsed  positions.  It  is  this 

distribution  of  orientations  about  the  common  bond  linking 

tetrahedra  that  has  a  profound  effect  on  the  third  peak.  Polk's 

model  also  has  a  density  within  2%  of  the  crystalline  density 

which  is  consistent  with  the  more  recent  experimental  r.d.f.'s 

Amorphous  germanium  samples  were  prepared  from  bulk  ingot 

by  resistance  heating  of  a  tungsten  boat  and  evaporation  onto 

-7 

1  mil  A1  foil.  The  vacuum  in  the  system  was  about  5  x  10  torr 
before  deposition  and  rose  to  about  10  ^  because  of  radiant 
heating.  The  thickness  of  the  foils  was  approximately  one  micron; 
six  layers  were  used  to  obtain  the  EXAFS  spectra.  These  samples 
were  studied  by  x-ray  dif fraction ,  electron  microscopy  and 
diffraction,  and  differential  scanning  calorimetry.  The  struc¬ 
tural  investigations  showed  diffraction  patterns  with  two  diffuse 

/  _1 
rings  at  %sine/x«  1.8l  and  3.22ft  .  Application  of  the 

Scherrer  fromula  to  the  line  width  gave  an  apparent  domain  size 

of  approximately  40ft.  Microscope  investigations  revealed  a  very 

find  40ft  wavelet  structure  on  top  of  a  300-500ft  "pebble  grain" 

texture.  These  observations  of  graininess  in  amorphous  films  are 

consistent  with  the  observations  of  others  and  do  not  mean  that 

•jg 

the  .  aterial  is  miorocrystalline  .  in  the  microscope  crystal¬ 
lization  could  be  induced  with  the  electron  beam  and  hot  stage 
examination  proved  a  first  crystallization  step  at  150-200°C. 

This  was  confirmed  with  differential  scanning  calorimetry  which 
showed  a  slow  exothermic  reaction  beginning  about  100-150°C  and 
continuing  until  the  rapid  onset  of  crystallization  about  450°C. 
These  calorimetric  observations  agree  with  the  earlier  work  of 
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Chen  and  Turnbull  .  Further  work  of  this  nature  is  being  carried 
out  to  correlate  the  structural  features  observed  with  EXAFS 
with  features  observed  with  these  other  techniques. 

The  crystalline  samples  consisted  of  a  pure,  finely  ground 
powder  (5  9's  pure  and  -325  mesh)  which  was  cast  in  Duco  cement 
and  spread  into  thin  layers  on  1  mil  A1  foil  with  a  glass  slide. 
Several  layers  were  used  to  give  an  effective  sample  thickness  of 
about  6u. 

The  EXAFS  data  for  crystalline  and  amorphous  Ge  are  shown 
in  fig.  12.  The  crystalline  data  is  the  smoothed  curve  of  fig. 

51).  The  amorphous  Ge  is  the  smoothed  curve  obtained  from  10  data 

runs.  As  can  be  seen  in  the  figure,  the  crystalline  data  has 

considerably  more  detail  than  the  amorphous  curve;  however,  the 
major  extrema  of  both  curves  coincide.  The  amorphous  curve  a 
appears  to  be  primarily  from  a  single  coordination  shell.  However, 
contributions  from  higher  shells  can  be  seen  above  600  eV. 

The  radial  structure  functions  (r.s.f.)  of  crystalline  and 
amorphous  Ge  obtained  from  the  Fourier  transform  of  the  data  are 
shown  in  fig.  13.  The  first  peak  on  both  curves  has  been  normal¬ 
ized  to  unity  since  only  the  relative  heights  of  the  peaks  are 

compared. 

In  crystalline  Ge  the  first  6  coordination  shells  are  clearly 
resolved.  The  expected  positions  of  the  coordination  shells  are 
given  in  Table  I  in  Chapter  2,  using  a  lattice  parameter  for  Ge 
of  5.65A,  while  the  noted  distance  in  the  figure  are  good  to  1% 
or  better.  It  should  be  emphasized  that  EXAFS  looks  at  crystal 
structure  differently  from  classical  crystallography  since  the 
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Figure  13 


Radial  structure  functions  of  crystalline  and 
amorphous  Oe  obtained  from  the  curves  cf  Fip;.  11. 

The  numbers  over  the  peak  indicate  the  measured  distance 
in  ft.  The  vertical  error  bar  (±lo)  in 
Fip;.  12,  1*1,  15,  16  was  estimated  by  replicate  analysis 
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structure  is  probed  radially  with  a  given  type  of  absorbing  atom 
at  the  origin.  The  transform  used  in  obtaining  fig.  13  (crystal¬ 
line  Ge)  was  calculated  without  the  phase  shift,  so  as  mentioned 
in  Chapter  3>  the  peaks  were  all  shifted  by  .28ft  to  the  expected 
position  given  in  Table  II.  Note  that  this  shift  is  in  good 
agreement  with  the  value  of  .27ft  obtained  from  the  linear  fit  of 
the  WKB  shift  in  Chapter  2. 

Another  parameter  that  was  obtained  from  the  crystalline 


data  is  y  .  As  mentioned  in  Chapter  3  if  N  ,  r  ,  and  o  are.  known 

J  J  J 

then  y  may  be  obtained  using  the  observed  A^'s  (noramlized  to  A^) 
and  eq.  3-5.  The  plot  obtained  from  the  data  is  shown  in  fig.  14. 
A  good  straight  line  is  obtained,  the  slope  of  which  is  y  =  .31ft"1 
which  corresponds  to  a  mean  free  path  for  scattering  into  an 
eigenstate,  1/  y  «=  3.2ft.  Because  this  value  is  of  the  order  of 


the  interatomic  spacing  and  because  the  short  range  structure  in 
crystalline  and  amorphous  Ge  is  so  similar  as  shown  in  fig.  13, 
it  is  assumed  that  the:  same  value  of  y  applies  to  amorphous  Ge. 

As  mentioned  previously,  this  assumption  is  critical  since  there 
is  no  way  to  independently  determine  y  for  an  amorphous  material. 

In  the  amorphous  r.s.f.  the  positions  of  the  first  two  peaks 
agree  exactly  with  the  crystalline  distance  while  the  third  peak, 
is  unresolvable  above  the  background. 

These  results  differ  from  earlier  results  which  found  that 


the  interatomic  distances  were  expanded  ^3%  from  the  crystalline 
.  41 

value  ,  and  are  more  in  agreement  with  mc-,e  recent  studies  on 

37  3  q 

Ge  and  Si  *  .  This  EXAPS  technique  is  not  sensitive  to  macro- 
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figure  14 

Plot  of  eq.  (3-5)  using  the  measured  crystalline  he  amplitudes. 
The  slope  of  the  curve  gives  y=.31. 
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or  micro-voids;  however,  a  uniform  distribution  of  atomic  sized 
vacancies  would  contribute  to  the  measured  values  by  causing 
distortions  in  the  average  atomic  distances.  Micro-voids  or  larger 
vacancy  clusters  were  not  present  since  electron  diffraction  did 
not  show  appreciable  small  angle  scattering. 

Since  the  first  distance  corresponds  to  the  tetrahedral  bond 
distance  and  the  second  distance  to  the  edge  of  the  tetrahedron, 
this  result  is  compatible  with  a  continuous  random  network  of  Ge 
tetrahedra,  i.e.  the  basic  tetrahedral  unit  of  crystalline  Ge  is 
preserved  in  amorphous  Ge  so  that  the  ideal  covalent  bonding  of 
each  Ge  atom  is  satisfied.  The  tetrahedra  are  not  linked  in  the 
regular  staggered  configuration  which  leads  to  the  diamond 
structure,  and  the  bond  angle  is  distorted  about  its  regular 
tetrahedral  value.  This  destroys  the  long  range  crystalline  order. 
This  is  indicated  by  the  broadened  and  reduced  second  peak.  The 
first  peak  in  the  crystalline  and  amorphous  curves  are  very  similar 
in  size  which  is  in  agreement  with  the  findings  of  Moss  and 
Gracezyk37  who  found  that  the  width  of  the  first  amorphous  peak 
in  Si  was  the  same  as  the  crystalline  width  and  thermal  in  origin. 
If  it  is  assumed  that  the  number  of  atoms  in  the  second  shell  and 
the  thermal  factor  are  the  same  as  in  crystalline  Ge,  then  by 
using  eq.(3-3)the  amplitude  difference  may  be  related  to  a  dif¬ 
ference  in  due  to  disorder  broadening.  By  eq.^3-7jit  is  found 

od2(A)  %  6  0TOT(C)  (5-D 

2  2 
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which  gives  a  disorder  smearing  of  .15ft  meaning  that  the  tetra¬ 
hedral  bonds  are  distorted  ±5°. 

/ 

D.  (let) 

fieO  is  isostructural  with  « -SiO^  and,  like  he ,  has  a  basic 
tetrahedral  unit  where  every  he  is  tetrahedrally  bonded  to  4 
oxygens  and  the  oxygen  atoms  are  shared  by  adjacent  tetrahedra. 
Previous  fine  structure  studies  have  shown  that  the  EXAWS  of 
amorphous  (Je02  is  very  similar  to  the  hexagonal  or  a  -quartz  form 
of  heO  .^2 >^3  This  crystal  structure  is  a  loose  stacking  of 
tetrahedra  in  a  hexagonal  unit  cell  with  3  tetrahedra  per  unit 
cell. 11 11  The  results  of  earlier  studies  on  the  structure  of 
amorphous  HeO  1,2  ’ 46  and  SiO^7  were  consistent  with  a  random 

network  of  OeO  tetrahedra.  Mozzi  and  Warren47  did  a  thorough 
study  of  the  radial  distribution  function  of  amorphous  Si02 
using  pair  functions  and  found  good  agreement  between  the  data 
and  a  continuous  random  network  model  of  Evans  and  King  . 

The  crystalline  HeO^-hex  with  the  a— nuartz  structure  was 
obtained  from  Fairmount  Chemicals.  Absorbers  were  m-de  by  grind¬ 
ing  the  powder,  casting  it  in  Duco  cement,  and  spreading  thin 
layers  on  1  mil  A1  foil.  The  amorphous  OeO^-a  was  made  bv  melting 

CeO  in  a  quartz  tube  and  quenching  into  water.  The  lack  of 
2 

crystallinity  was  verified  by  x-ray  diffraction  which  showed  only 
a  series  of  broad,  smooth  maxima.  The  amorphous  absorbers  were 
prepared  in  the  same  fashion  as  the  crystalline  Ge02,  and  the 
thickness  of  both  absorbers  was  ^10  u. 

A  comparison  of  crystalline  and  amorphous  ReOg  is  shown 


in  fig.  15. 
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Figure  15 


Comparison  of  crystalline  and  amorphous  CeO^  EXA RC 


The  transforms  of  the  experimental  EXAFS  data  is  shown  in 

fig.  1 6.  This  is  data  for  only  the  Ge  edge  since  the  oxygen 

edge  is  beyond  current  instrumental  capability.  Both  curves 

have  the  first  peak  normalized  to  unity.  When  the  transforms 

were  taken  peaks  near  r  =  0  were  observed  which,  depended  on  where 

the  data  was  terminated  near  k  =  0.  The  data  also  contained 

termination  ripple  of  periodicity  as  discussed  in  Chapter  3. 

Typically  several  transforms  were  done  with  various  ranges  of 

data  in  order  to  get  transforms  which  were  nearly  flat  near 

r  =  0  and  to  identify  termination  ripple.  A  typical  transform 

curve  for  GeO  -hex  is  shown  in  fig.  17  with  residual  structure 
2 

and  termination  ripple  left  in,  and  the  peaks  shifted  by  "a” 
from  their  crystalline  distances.  The  final  curves  (fig.  lo) 
have  had  the  residual  structure  near  r  =  0  and  termination  ripple 
graphically  removed.  Work  is  currently  underway  to  improve  the 
data  analysis  techniques  to  eliminate  the  structure  added  by 
the  transformation  methods.  To  compare  the  crystalline  data  with 
the  expected  positions  and  amplitudes  of  the  peaks  from  knowledge 
of  the  crystal  structure,  a  program  was  developed  which  can 
calculate  the  number,  types,  and  distances  of  atoms  from  a  given 
atom  placed  at  the  origin.  This  program  is  explained  in  detail 
in  Appendix  B,  and  the  resulting  structural  information  is  listed 
in  Table  III.  The  lines  in  fig.  16  show  the  expected  peak 
positions  and  their  amplitudes  calculated  from  eq.^3-3^  In  the 
expression  for  the  amplitudes  shell  to  shell  variations  of  a 

J 

were  ignored  and  the  difference  in  scattering  factors  was  accounted 
for  by  assuming  that  each  factor  was  proportional  to  the  number 


SeGe-C 
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RAW  TRANSFORM  OF  GeO.-hex 


Figure  17 


Transform  of  GeO^-hex  data  showing 


termination  rinole  and  structure  near  r 
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Table  III 

EXAPS  Structural  Data  for  QeO^,,  GeSe ,  SeGe 


GeO  GeSe  SeGe 

2 


N 

J 

r  (ft) 

J 

atom  type 

N 

J 

r  (ft) 

J 

atom  type 

N 

J 

r  (ft) 

.1 

atom  type 

2 

1.63 

0 

1 

2 . 56 

Se 

1 

2.56 

Ge 

2 

1.65 

0 

2 

2.59 

Se 

2 

2.59 

Ge 

4 

3.13 

Ge 

2 

3.31 

Se 

2 

3.31 

Ge 

2 

3.56 

0 

1 

3.35 

Se 

1 

3.35 

Ge 

2 

3.62 

0 

2 

3.39 

Ge 

4 

3.68 

Se 

2 

3.70 

0 

2 

3.04 

Se 

2 

3.69 

Se  ' 

2 

*1.01 

0 

4 

3.94 

Ge 

2 

3.84 

Se 

2 

4.08 

0 

4 

4.04 

Ge 

1 

3.94 

Ge 

2 

4.17 

0 

2 

4.39 

Se 

2 

4.39 

Se 

2 

4.24 

0 

2 

4.62 

Se 

2 

4.62 

Ge 

2 

4.26 

0 

1 

4.65 

Se 

1 

4.65 

Ge 

2 

4.31 

0 

2 

4.72 

Ge 

2 

5.10 

Ge 

2 

4.42 

Ge 

2 

5.10 

Se 

3 

5.56 

Ge 

2 

4.52 

Ge 

3 

5.51 

Se 

2 

5.51 

Ge 

2 

4.58 

0 

2 

5.54 

Se 

4 

5.74 

Se 

2 

4.64 

0 

4 

5.84 

Ge 

8 

5.84 

Se 

2 

4.68 

0 

2 

6.01 

Se 

2 

6.01 

Ge 

8 

4.97 

Ge 

2 

6.02 

Se 

2 

6.02 

(i6 

5.19 

0 

2 

6.03 

Se 

2 

6.03 

Ge 

2 

5.29 

0 

4 

6.28 

Ge 

2 

6.37 

Ge 

2 

5.49 

Ge 

2 

6.37 

Se 

2 

6.41 

Ge 

2 

5.64 

Ge 

4 

6.41 

Se 

6 

6.58 

Se 

4 

5.74 

0 

2 

6.71 

Se 

2 

6.71 

Ge 

2 

5.95 

Ge 

2 

6.73 

Se 

2 

6.73 

Ge 

2 

5.97 

0 

5 

6.75 

lSe/4Ge 

1 

6.75 

Ge 

2 

5.99 

0 

4 

6.78 

Ge 

8 

6.99 

Se 

2 

6.02 

0 

2 

7.16 

Se 

2 

7.16 

Ge 

2 

6.05 

0 

2 

7.20 

Ge 

4 

7.23 

Se 

4 

6.08 

0 

2 

7.34 

Se 

2 

7.34 

Ge 

2 

6.12 

0 

2 

7.35 

Se 

2 

7.35 

Ge 
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of  electrons  as  shown  In  ChaDter  2,  so  that  f  =  *lf  .  To 

he  0 

obtain  a  good  agreement  between  the  expected  crystalline  peak 
positions  and  the  observed  peak  positions  all  of  the  peaks  in 
the  r.s.f.  had  to  be  shifted  by  +.28?.  which  is  the  same  phase 
shift  factor  found  in  crystalline  he. 

Comparing  the  crystalline  and  amorphous  data,  the  first 

peak,  which  is  the  tetrahedral  he-0  bond,  moved  in  slightly 

from  1.65?  to  l.oO?.  The  second  peak  which  is  the  he-he  distance 

in  adjacent  tetrahedra  decreased  from  3.15?  to  3.07?.  Although 

it  is  possible  that  some  of  this  difference  might  be  due  to 

differences  in  the  phase  shift  between  heO  -hex  and  heO  -  a 

2  2  5 

contrary  to  our  assumption  of  equal  phase  shifts  in  the  materials 
this  was  not  considered  as  likely  as  a  change  in  bond  length 
since  the  k  dependence  of  the  Ge  phase  shifts  in  crystalline 
Ge ,  OeO^  and  GeSe  was  found  to  be  .28^.03?  independent  of  the 
crystal  structure  or  whether  the  bond  is  covalent  or  ionic. 

'  h.is  means  that  the  Ge— 0— Ge  bond  angle  changes  from  its  value 
of  1^°  in  the  crystalline  case  to  liJ7°±5°  in  the  amorphous  case. 
The  25£  decrease  in  peak  amplitude  can  be  attributed  to  a  pro¬ 
portionate  decrease  in  the  number  of  atoms  or  disorder  smearing, 
or  both.  If  amorphous  GeO^  is  a  continuous  random  network  of 
tetrahedra,  then  the  decrease  is  more  likely  due  to  disorder 
smearing  as  in  Ge.  The  third  GeO^-a  peak  at  3.3^?  is  signifi  - 
cantlv  different  from  GeO^-hex  and  from  the  amplitude  of  the 
peak  it  could  be  due  to  1  Ge  atom  or  4  0  atoms  if  there  is  no 
disorder  smearing.  Since  this  distance  is  greater  than  that 
possible  for  a  Ge-Ge  distance  in  adjacent  oxygen  sharing  tetra- 
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hedra,  it  is  tentatively  Identified  with  the  median  oxygen 

atoms  in  adjacent  tetrahedra  and  suggests  that  the  alignment  of 

oxygen  atoms  is  more  symmetric  (like  the  S-cristobalite 

48 

structure  )  in  GeOg-a  than  in  GeO^-hex.  These  results  suggest 

that  GeO  is  a  continuous  random  network  of  tetrahedra  which  is 
2 

more  ordered  than  for  Ge ,  and  are  in  general  agreement  with  the 
earlier  studies  although  these  results  surpass  others  in  reso¬ 
lution  and  detail. 

C .  GeSe 

The  final  system  studies  was  GeSe  which  was  chosen  to 
illustrate  the  power  of  the  EXAPS  technique  to  probe  the  structure 
around  individual  atomic  species.  The  crystal  structure  of  GeSe 
is  highly  distorted  NaCl  structure  of  the  SnS  type'^’^1  with  each 
atom  having  6  neighbors  of  the  opposite  kind  (3  at  ^2. 58ft  and  3 
at  'v3«33ft)«  However,  beyond  the  second  neighbors  the  structure 
about  a  Ge  atom  and  a  Se  atom  differ.  Unlike  Ge  and  GeO^,  there 
is  no  easily  identifable  basic  structural  unit  in  either  the 
crystalline  or  amorphous  case. 

Crystalline  GeSe  was  prepared  by  heating  equiatomic  portions 
of  pure  Ge  and  Se  in  a  solid  quartz  tube  for  several  hours  at 
900°C  and  quenching  in  water.  The  crystal  structure  was  verified 
from  the  x-ray  diffraction  pattern  which  agreed  with  the  listed 
pattern  in  the  ASTM  files.  Absorbers  were  prepared  by  grinding 
the  GeSe  into  a  fine  powder  and  mixing  it  with  Duco  cement  and 
spreading  thin  layers  onto  A1  foil.  Amorphous  samples  were 
prepared  by  vacuum  evaporation  under  the  same  conditions  as 
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amorphous  Ge.  The  problem  of  non-stoichiometric  sublimation  of 
GeSe  was  recognized  }  and  since  a  multiple  furnace  technique 
was  not  available,  various  sized  pieces  of  GeSe  were  sublimed 
simultaneously  in  order  to  get  a  uniform  deposition.  Measure¬ 
ments  of  the  relative  x-ray  absorption  coefficients  proved  that 
the  resultant  films  were  within  3$  of  the  composition  of  the 
initial  evaporants  although  this  does  not  prove  that  the  material 
is  homogeneous.  Some  differential  scanning  calorimetry  was  done 
which  showed  a  single  exothermic  crystallization  cealc  at  about 
350° C  which  indicated  that  the  material  is  homogeneous  and  not 
phase  separated  Ge  or  Se,  which  crystallize  at  450°C  and  70°C 
respectively.  The  thickness  of  both  the  crystalline  and  amorphous 
GeSe  absorbers  was  'vlOp. 

The  K  edge  absorption  fine  structure  for  both  the  Ge  and  Se 
atoms  in  amorphous  GeSe  was  shown  in  fig.  1  and  a  comparison  of 
the  crystalline  and  amorphous  EXAFS  for  both  edges  is  shown  in 
fig.  18 .  The  differences  in  the  EXAFS  can  be  seen  although  a 
close  inspection  would  show  that  the  basic  periodicity  of  the 
major  extrema  is  the  same.  The  radial  structure  functions  of 
both  crystalline  and  amorphous  GeSe  is  shown  in  fig.  19  a  and  b. 
Note  that  the  r.s.f.’s  are  relative  to  a  Ge  and  an  Se  atom, 
respectively.  As  in  fig.  16  the  lines  indicate  the  positions  of 
the  atoms  in  the  crystal  relative  to  either  a  Ge  or  Se  atom  cal¬ 
culated  from  the  computer  program  of  Appendix  C  with  the  amplitudes 

from  eq.  3-3  ignoring  possible  differences  in  the  o  's.  The 

J 

crystalline  data  shows  that  the  experimental  r.s.f.'s  from  each 
edge  agree  quite  well  with  the  predicted  structure  and  that  the 


SeGe-A 
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Xr/  =  d/oi)  u* 


Comparison  of  crystalline  and  amorphous 


Figure  19 


Radial  structure  functions  for  crystalline  and  amornhous  fleSe 
with  respect  to  a  Oe  atom  (l?a)  and  an  Se  atom  (19b)  at  the  origin 
Phe  lines  indicate  the  positions  and  amplitudes  of  the  crystalline 
peaks  expected  about  each  kind  of  atom  from  a  knowledge 
of  the  crystal  structure  and  eq.  (3-3). 
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Comparison  of  the  rsf's  about  the 
He  and  Se  atom  in  amornhous  CeSe. 
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at  2.35ft  which  gives  1  atom  at  2.80ft  using  eq.  3-3.  Since  the 
ideal  GeSe  covalent  bond  is  2.37ft,  it  appears  that  in  going 
from  the  crystalline  to  the  amorphous  material  the  material  has 
changed  from  a  more  to  a  less  distorted  covalent  bonding  arrange¬ 
ment  for  the  nearest  neighbor. 

In  fig.  20  the  splitting  around  the  Se  atom  is  more  clearly 

resolved  and  may  illustrate  the  tendency  of  Se  to  link  two  nearby 
Ge  atoms  into  a  chain  as  the  other  nearby  Ge  atom  is  forced  away. 
The  less  resolved  splitting  about  the  Ge  atom  may  be  due  to  bond¬ 
ing  to  another  Ge  atom  (perhaps  at  the  Ge-Ge  bond  distance  of 
2.A5ft)  which  would  indicate  the  tendency  of  Ge  to  form  tetra- 
hedrally  directed  bonds.  A  satisfactory  structural  unit  incor- 
oorating  these  features  and  satisfying  all  bonding  requirements 
has  not  yet  been  visualized. 


Chapter  6 

Discussion  and  Conclusions 


Although  EXAFS  has  been  shov/n  to  arise  from  a  complex  Inter¬ 
action  of  the  ejected  photoelectron  with  the  surrounding  atoms,  it 
is  measured  by  simnly  determining  the  probability  of  an  x-ray 
photon  being  absorbed  at  a  given  energy.  This  probability  is 
related  to  the  overlap  of  the  initial  K  electron  wave  function 
with  the  final  ejected  photoelectron  wave  function.  If  the  final 
photoelectron  wave  function  contains  significant  scattering 
contribution  from  atoms  nearby  to  the  absorbing  atom,  then  oscil¬ 
lations  in  the  absorption  coefficient  as  a  function  of  energy 
will  appear.  Thus,  oscillations  in  the  absorption  probability 
(which  are  measured  as  variations  in  1/1°  as  a  funotion  of  energy) 
are  determined  by  the  location  of  the  surrounding  atoms.  This 
interaction  is  "sensed"  by  the  x-ray  photoelectron  during  the 
time  of  its  interaction  with  the  K  shell  electron.  The  quantum 
mechanical  calculation  is  straight  forward  once  the  physical 
mechanism  is  understood.  Treatment  of  the  surrounding  atoms  as 
a  system  of  weak  point  scatterers  simplifies  the  calculation  of 
the  scattered  waves  yet  retains  the  physics  of  the  problem.  The 
final  result  (eq.  2-33)  contains  the  energy  dependence  of  the 
scattering,  a  temperature  factor,  and  a  range  term  which  allows 
calculation  of  spectra  in  quantitative  agreement  with  the  observed 
EXAFS  for  three  different  crystal  structures.  Eq.  2-33  is  suf¬ 
ficiently  general  that  it  can  be  applied  to  any  system  where  EXAM’S 
is  observed  (crystalline  and  non-crystalline  solids,  liquids,  and 
polyatomic  gases)  and  for  either  unpolarized  or  polarized  x-rays. 
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For  unpolarized  x-rays  (e  ^  =  e  ^  a  e  or  for  materials  with 

^  y  ^ 

cubic  symmetry  (x2=y2=z2),  the  result  simplifies  to  eq.  2-3*1. 

i  i  i 


However,  this  result  may  also  be  used  to  explain  the  orientation 

dependence  of  EXAFS  observed  with  polarized  x-rays  and  single 

crystals  of  non-cubic  materials  in  a  more  direct  way  than  has 

been  done  previously.  Future  work  is  planned  to  use  the  point 

scattering  model  to  calculate  the  L  and  L  spectra  which  are 
b  II  III 


dominated  by  the  2p  electrons  making  transitions  to  d  symmetric 
states  in  the  continuum.  These  calculations  are  complicated  by  a 
more  complex  angular  dependence  for  both  the  initial  and  final 
electron  wave  function. 

The  most  important  result  of  the  point  scattering  theory 
shows  that  EXAFS  can  be  Fourier  transformed  to  give  a  radial 
structure  function  for  each  kind  of  atom  in  a  material.  The 
r.s.f.  contains  information  on  the  number  and  kinds  of  atoms  and 
their  distribution  about  the  coordination  distances.  The  r.s.f. 
is  superior  to  the  radial  distribution  function  in  several  ways. 
The  main  advantage  is  the  radial  distribution  function  only 
calculates  a  two-body  correlation  function,  whereas  the  r.s.f. 
determines  the  average  structure  about  each  different  kind  of  atom 
in  a  multicomponent  system.  This  information  (note  the  GeSe  data) 
can  be  essential  for  understanding  the  structure.  Another 
advantage  is  that  the  r.s.f.  has  potentially  twice  the  resolving 
power  of  the  r.d.f.  since  there  are  twice  as  many  oscillations  in 
the  EXAFS  data  over  the  same  range  of  k  space.  And,  finally,  the 
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natural  definition  of  the  r.s.f.  as  a  sum  of  Oaussian  curves 
from  each  coordination  shell  allows  the  height  and  width  para¬ 
meters  to  be  measure^  for  each  shell,  while  the  r.d.f.  [defined 
as  4irr2p(r)]  determines  the  positions  of  the  atomic  shells  as 
fluctuations  about  the  average  (parabolic)  density  curve.  Fitting 
Haussian  curves  to  determine  the  area  of  the  r.d.f.  neaks  is 
usually  only  possible  for  the  first  peak  and  even  then  may  be 
ambiguous  because  of  contributions  from  different  kinds  of  atoms 
in  a  multicomponent  system. 

Many  of  the  disadvantages  of  the  EXAFS  technique  can  be 
overcome  by  further  developments.  These  include  improvements  in 
data  analysis  to  lessen  termination  effects  in  the  data,  develop¬ 
ment  of  a  theory  for  the  L^and  L  edges  to  obtain  structural 
information  from  those  materials  where  the  K  shell  is  too  high 
in  energy,  acquisition  of  a  vacuum  double  crystal  x-ray  spectro¬ 
meter  to  extend  the  range  of  material  whose  K  or  L  edges  may  be 
measured  (this  includes  the  whole  periodic  table),  and  improved 
data  acquisition  techniques  to  speed  up  and  simplify  the  data 
Gathering  and  analysis.  ^ 

There  remain  some  inherent  disadvantages  of  the  EXAFS 

method:  The  amorphous  data  must  be  compared  to  the  crystalline 

data  in  order  to  account  for  parameters  such  as  e  ,  y ,  n ,  and  the 

T 

amplitude  normalization  for  the  material.  A  more  serious  limitation 
is  that  in  multicomoonent  samples  the  relative  absorption  of 
each  element  whose  EXAFS  is  to  be  studies  must  be  large  enough 
so  that  the  fine  structure  can  be  resolved,  (e.g. — a  few  atomic  % 
of  a  light  element  in  a  heavy,  element  matrix  could  not  be  measured.) 
This  is  the  main  experimental  limitation,  but  it  still  leaves  a 
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wide  range  of  systems  which  can  be  studied. 

Some  tentative  conclusions  on  amorphous  structures,  can  be 
drawn  from  the  three  systems  studied  here.  Two  of  the  systems 
(fie  and  GeO  )  were  examples  of  a  continuous  network  of  tetrahedra 
In  both  cases  the  tetrahedral  unit  which  is  the  basis  of  the 
crystal  structure  is  preserved  in  the  amorphous  stystem.  In  Ge 
each  atom  forms  an  ideal  covalent  bond  with  4  neighbors  at  a 
distance  of  2.45ft,  while  in  Ge02  each  Ge  has  4  tetrahedrall.v 
bonded  oxygens  with  each  oxygen  acting  as  a  bridge  between  tetra¬ 
hedra  by  bonding  with  two  Ge.  The  Ge-0  bond  is  partially  ionic 
with  a  length  of  1.60ft.  The  networks  built  up  from  these  tetra¬ 
hedra  are  substantially  different.  In  Ge  adjacent  tetrahedra 
are  oriented  randomly  about  the  common  Ge-Ge  bond  which  leads 
to  a  very  rapid  reduction  in  order  beyond  the  second  distance. 

The  GeO  network  differs  from  Ge  since  the  oxygen  atoms  must  be 
shared  by  adjacent  tetrahedra,  so  that  the  order  in  this  system 
extends  out  beyond  the  neighboring  tetrahedra.  In  GeSe  the 
situation  is  much  different.  The  crystalline  state  is  a  highly 
distorted  NaCl  structure  with  each  atom  having  3  neighbors  of 
the  opposite  kind  at  -v2.58ft  and  3  neighbors  of  the  opposite  kind 
at  3.33R.  These  bonds  are  highly  distorted  from  the  ideal  co¬ 
valent  or  ionic  bond.  In  the  amorphous  material  the  structure 
changes  so  that  2  of  the  atoms  move  into  a  distance  of  2.35ft 
(which  is  the  GeSe  convalent  bond  length),  while  the  other  near¬ 
neighbor  atom  moves  out  to  2.80ft.  Around  the  Ge  atom  there  is 
also  evidence  of  a  4th  atom  at  2.45ft  (the  Ge-Ge  distance).  All 
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of  these  results  show  that  each  atom  tries  to  satisfy  its  ideal 
bonding  conditions  in  the  amorphous  state  (i.e.  2  atoms  in  a 
chain  in  Se,  and  4  tetrahedrall.v  bonded  neighbors  in  Ge).  We 
make  the  following  general  postulate  for  amorphous  structures: 

In  covalently  or  ionically  bonded  amorphous  systems  when  the 
requirements  of  long  range  order  are  removed,  the  near  neighbor 
structure  is  determined  when  each  atom  achieves  its  ideal  bonding 
state  to  the  extent  possible  considering  the  composition  of  the 
material.  This  postulate  supports  the  earlier  models  of  satis¬ 
faction  of  each  atom's  covalent  bonding  requirement  in  amorphous 
semiconductors  by  Mott5'1  and  Cohen>  Prltzsche>  and  Ovshlnsky . 55 

Further  work  is  planned  on  certain  compositions  in  the  Ge-As-Se 
system  to  verify  this  statement. 

In  summary,  the  main  conclusion  of  this  work  is  that  K  edge 
EXAFS  is  given  by 

xoo  =  -kf (k)  y 


N,e“Yrj  -2o,2k2 


J  sin  (2kr  +  2n(k))  ^  ^ 


(e  2x, 


+  ev  yi 


ezV> 


which  not  only  allows  calculation  of  fine  structure  spectra  which 
are  in  good  agreement  with  observed  spectra  for  any  x-ray  polar¬ 
isation  system  where  the  structure  is  known,  but  whose  Fourier 
transform  allows  structure  determination  in  systems  where  the 
structure  is  unknown.  Application  to  several  systems  has  clearly 
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demonstrated  that  it  is  a  powerful  technique  with  some  unique 
capabilities  with  respect  to  the  standard  scattering  techniques. 

-14  .  . 

The  EXAFS  technique  shows  promise  of  becoming  the  primary  method 
for  structure  determination  in  amorphous  materials. 


Appendix  A 
Filtering  Methods 


The  absorption  spectrum  obtained  from  the  x-ray  data  contains 
three  components;  (1)  the  K  absorption  edge  followed  by  a  mono- 
tonically  decreasing  absorption  coefficient  above  the  edge;  (2) 
the  EXAFS  added  to  this  above  the  edge,  and  (3)  the  random  statis¬ 
tical  counting  noise  superimposed  on  all  of  the  data.  The  EXAFS 
must  be  separated  from  the  other  two  components .  This  could  be 
done  by  fitting  some  functional  form  to  the  usual  photoelectric 
spectrum,  subtracting  it  from  the  data,  and  then  using  a  running 
average  of  data  points  to  reduce  the  effect  of  the  noise.  Ihe 
major  difficulty  is  that  the  functional  form  of  each  edge  is  not 
well  known  and  would  require  an  extensive  fitting  effort  for  each 
absorption  edge. 

The  approach  used  here  entails  filtering  in  the  large  in  which 
the  data  is  Fourier  transformed  and  the  frequency  components  of 
the  desired  signal  are  extracted  from. the  whole  spectrum  by  filter¬ 
ing  out  the  undesired  components  and  then  retransforming  the  filter¬ 
ed  data.  This  method  is  particularly  suited  for  the  EXAFS  tech¬ 
nique  since  the  first  component  contains  only  very  low  frequency 
components,  the  EXAFS  signal  (component  2)  is  limited  to  a  small 
band  of  frequencies  separated  from  the  low  frequency  end,  while 
the  random  noise  (component  3)  bas  contributions  over  the  entire 
domain.  The  photoelectric  background  is  removed  by  first  subtract¬ 
ing  the  lowest  frequency  components  from  the  transformed  data. 
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Although  some  small  contribution  from  higher  frequencies  may  still 
exist  In  the  data,  close  Inspection  of  the  EXAPS  obtained  using 
this  method  showed  no  residual  effects  of  the  background. 

Since  the  EXAPS  components  are  limited  to  the  middle  range  of 
the  total  domain, the  noise  Is  removed  by  filtering  out  all  higher 
components.  This  is  shown  In  fig.  21  where  the  arrows  indicate  the 
high  and  low  cuts.  Although  the  smoothed  EXAPS  contains  a  small 
contribution  from  the  components  of  the  noise  within  the  same  band 
as  the  signal,  this  is  only  important  at  the  upper  end  of  the 
spectrum  as  can  be  seen  in  fig.  5b.  When  the  transformation  tech¬ 
niques  of  Chapter  3  are  applied  to  the  data  the  noise  is  left  in 
to  assure  that  all  of  the  signal  is  also  included. 

The  actual  computer  algorithm  used  to  transform  the  data  is 
the  HARM  program  from  the  IBM  Scientific  Subroutine  Package  which 
is  based  on  the  Fast  Pourier  Transform  method  originally  developed 
by  Cooley  and  Tukey.56  This  program  requires  that  the  data  be  taken 
in  equally  spaced  intervals  in  one  domain  and  that  the  number  of 
points-be  a  power  of  2  (usually  210  =  1024  points;  zeros  can  be 
appended  to  get  the  proper  number  of  points).  The  data  is  Pourier 
analyzed  in  terms  of  the  real  and  imaginary  components  of  512 
frequencies.  The  lowest  frequency  has  a  wavelength  of  twice  the 
range  of  the  data  and  the  remaining  frequencies  are  the  next  511 
harmonics  of  this  fundamental.  The  filtering  is  then  done  in  terms 
of  .this  finite  set  of  frequencies. 

lhe  actual  data  is  taken  in  equal  steps  in  angle  ( 2  ©) .  The 
'FT  filtering  to  remove  the  photoelectron  background  is  done  with 
this  data.  However,  during  the  transformation  to  obtain  structural 
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Figure  21 

i'-xamnle  of  a  frequency  spectrum  of  the  raw  data  using  the 
fast  Fourier  transform.  The  arrows  show  the  high  and  low 
frequency  cutoffs.  Only  70  of  the  512  components  are  shown; 
all  higher  components  are  negligible. 
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information,  an  interpolation  routine  is  used  to  generate  a  set 
of  data  which  is  equally  spaced  in  k  and  then  Fourier  transformed, 
as  in  Chapter  3,  to  obtain  the  radial  structure  function. 


Appendix  B 

Atomic  Structures  for  EXAPS 

In  the  point-scattering  theory  of  EXAPS  the  positions  of  the 
atoms  are  located  by  the  ejected  photoelectron  as  it  propagates 
radially  outward  from  the  absorbing  atom.  In  crystalline  materials 
the  surrounding  atoms  are  located  from  a  given  atom  at  a  discrete 
series  of  distances  known  as  coordination  shells.  This  is  a 
different  way  of  looking  at  crystal  structure  than  in  classical 
crystallography  which  gives  the  structure  in  terms  of  a  unit  cell 
which  may  or  may  not  have  a  particular  atom  at  the  origin,  and,  in 
fact,  the  crystal  structure  in  terms  of  coordination  shells  is 
only  listed  for  a  few  of  the  simplest  crystal  structures  (simnle 
cubic,  fee,  bcc,  dia  cubic,  hep). 

To  calculate  EXAPS  from  eq.  2-3*1  for  a  known  crystal  structure 
or  to  determine  the  location  and  amplitudes  of  the  peaks  in  a 
crystalline  r.s.f.  using  eq.  3-3,  a  general  computer  program  was 
developed  which  can  compute  the  structure  in  terms  of  coordination 
shells  about  a  given  atom  if  the  unit  cell  is  known.  This  program 
will  work  for  any  symmetry,  and  for  any  number  of  atoms  in  the  unit 
cell  and  is  limited  only  by  the  oore  size  of  the  computer  used. 

The  input  required  is  (1)  the  lengths  and  angles  between  the 
coordinate  axes  in  the  unit  cell  (a,b,c:  c(,  g,  y  in  the  usual  crystal¬ 
lographic  notation),  (2)  the  number  of  kinds  of  atoms,  (3)  the 
number  of  atoms  of  each  kind  in  the  unit  cell,  (4)  the  maximum 
number  of  lattice  translations  along  each  coordinate  axes,  and  (5) 
if  there  is  no  atom  at  the  origin,  the  displacement  needed  to  bring 
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a  particular  atom  (corresponding  to  the  absorbing  atom)  to  the 
origin.  Because  the  program  tallies  each  atomic  position,  including 
the  positions  of  the  translations  of  the  unit  cell,  core  limitations 
permit  only  those  atom  positions  in  a  volume  within  about  10  near 
neighbor  distances  in  each  direction  to  be  determined.  This  is  not 
a  severe  restriction  since  EXSPS  only  probes  about  10*  from  the 
absorbing,  atom,  so  that  in  atoms  with  large  unit  celis  (e.g.  with 
72  atoms  in  a  unit  cell)  only  one  or  two  translations  in  each 
direction  are  needed  to  find  all  of  the  distances  of  interest  for 

EXAM’S. 

The  computer  program  calculates  the  x,y,  and  z  coordinates 
for  all  of  the  atoms  of  a  given  type  relative  to  a  particular  atom 
at  the  origin  using  a  nest  of  DO  loops.  The  distances  of  the 
atoms  from  the  origin  are  then  computed  using  the  general  formula 

r  _  (x  2  +  y  2  +  2  2  +  2xy  cos  Y  +  2yz  cos  a  +  2xz  cos  3)  (B-l) 

where  a  ,  '3  ,  and  y  are  the  angles  between  the  coordinate  axes.  The 
program  prints  out  the  type  of  atom,  which  atom  of  that  type  it  was 
the  translations  made  along  each  coordinate  axes  and  the  total 
distance  (using  eq.  B-l),  and  a  histogram  showing  the  distances  and 
number  of  atoms  at  the  distance.  If  the  compound  of  interest  has 
two  types  of  atoms  (say  A  and  B  with  A  being  the  absorbing  atom), 
the  program  first  calculates  the  A-A  distances  and  histogram,  then 
the  A-B  distances  and  histogram,  and  finally  a  total  histogram. 

The  same  thing  could  then  be  done  for  the  B  atoms  by  moving  a  B 
atom  to  the  origin  of  the  coordinate  system. 
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This  program  is  also  useful  in  materials  where  there  may  be 
several  non-equivalent  positions  for  a  given  type  of  atom.  The. 
structure  for  each  separate  position  can  be  found  by  moving  each 
different  position  to  the  origin.  The  structure  used  to  calculate 
EXAFS  is  then  the  average  structure  from  the  different  positions 
of  the  same  atom  since  the  EXAFS  just  reflects  the  average 
structure  about  a  particular  kind  of  atom. 

Typical  data  from  this  program  has  been  shown  in  Table  III 
and  a  listing  of  the  program  is  given  in  Appendix  C. 


Appendix  C 

Computer  Program  Listings 
A  EXAFS  Atomic  Structure  Program 


C  EXAFS  STRUCTURE  PROGRAM  USING  EQ  (B-l) 
C 


C 

C 

C 


C 

C 

C 


REAL  DAZ(3,10) ,RHO(2 ,4 
REAL  UX(3,2H),UY(3,2H), 
DIMENSION  NOPOCK (1001) , 
DATA  IASTER/1H*/ 

DATA  IBLANK  /1H  / 

MAXSUB-501 

MAXSUB-901 


9,9,9) 

UZ(3,2i»),DAX(3,2i»),DAY(3,2i<) 

N0T0TP(1001),N(5) 


READ  INPUT  DATA 

READ(5,901)  A,3,C,AL,BE,GA 

READ( 5 ,902 )  MH,MK,ML,NAS,N(1),N(2) 

NA*N ( 1 ) NB=N ( 2 ) 

READ  ATOMIC  COORDINATES 


C 

C 

C 


DO  100  J-l.NA 

HFAD(5,9^)  UX(1,J)  ,UY(1,J),UZ(1,.T) 

100  CONTINUE 

IF(NB.EQ.O)  GO  TO  200 
DO  200  K=1,NB 

READ(5,90iJ)  UX(2,K)  ,  UY(2,K),UZ(2,K) 
200  CONTINUE  *  *  1  *  ' 


RAD=  3.1^15926/180. 

ALPIIA=AL#RAD 

BETA=BE*RAD 


GAMMA =GA#RAD 

MAH=2#MH+1 

MAK=2#MK+1 


MAL=2»ML+1 


INITIALIZE  HISTOGRAM 


DO  205  1=1,MAXSUB 
N0T0TP(I)*0 
205  CONTINUE 

DO  1  1=1, NAS 
NA-N(I) 

WRITE ( 6, 9 02) NAS, NA 
DO  210  IPOCK=l ,MAXSUB 
NOPOCK (I POCK )=0 


210  CONTINUE 

WRITE (6, 905) 

Q 

C  FORM  X  COORDINATES  FOR  ALL  ATOMS 
C 

DO  2  IH-1,MAH 

X=IH-MH-1 

ix=x 

DO  10  JA=1,NA 

10  DAX(I , JA)=(X+UX( I , JA) )*A 

Q 

C  FORM  Y  COORDINATES  FOR  ALL  ATOMS 
C 

DO  3  IK-l.MAK 

Y-IK-MK-1 

IY*Y 

DO  20  JB=1,NA 

20  DAY(I,JB)=(Y+UY(I,JB))«B 

Q 

C  FORM  Z  COORDINATES  FOR  ALL  ATOMS 
C 

DO  4  IL-l.MAL 

Z»IL-ML-1 

IZ=Z 

DO  30  JC-1,NA 
DAZ(I,JC)=( Z+UZ (L, JC) ) *C 
DX-DAX(I,JC) 

DY=DAY (I ,JC) 

DZ=*DAZ ( I ,  JC) 

C  TRANSLATE  COORDINATES  SO  THAT  A  PARTICULAR  ATOM  IS  AT  THE  ORIGIN 
C 

DX=DX-. 226*A 
DY=DY-.226#B 
DZ=DZ-.226*C 

Q 

C  GENERATE  DISTANCES  USING  EQ  (B-l) 


RHO(I  JC,IH,IK,IL)=SQRT(DX*DX+DY*DY+DZ*DZ+2.*DX*DY*COS(GAMMA) 

1  +2 . #DY*DZ*COS( ALPHA)+2 . #DX*DZ#COS (BETA ) ) 

DIST*RHO ( I , JC , IH , IK . IL) 

NSUB=IFIX(DIST*100.)+1 
MARK-I ASTER 

IF(NSUB.GT.O  .AND.  NSUB.LE ,MAXSUB)GO  TO  220 
MARK-I BLANK 
GO  TO  222 

220  NOPOCK  (NSUB)  =NOPOCK  (NSUB)+.l 
NOT OTP ( NSUB ) =NOTOTP ( NSUB ) +1 


3  WRITE  OUTPUT 

'222  WRITE(6  , 906)1, IX,IY,IZ,JC,DIST, MARK 
30  CONTINUE 


o  o  o 


105 


GENERATE  THE  HISTOGRAMS 

WRITE (6, 90 8) 

90 8  FORMAT  (1111) 

DO  230  IHIST-1, MAXSUB 
BOT-FLOATC IHIST-1 ) * . 01 
TOP= FLOAT ( IHIST ) * . 01 
NPHIST-NOPOCK ( IHIST) 

IFCNPHIST  .LE.  110  )GO  TO  225 
WRITE (6, 810) 30T, TOP, I 
GO  TO  230 
225  CONTINUE 

IFCNPHIST. GT. 0) WRITE ( 6, 809)3OT, TOP, NPHIST, I, (IASTER, IAST-1 , 
1NPIIIST) 

809  F0RMAT(IX,F5.2,F7.2,15,2X,I1,110A1) 

810  FORMAT ( IX, F5 . 2 ,F7 . 2  ,7X,I1 , '  TOO  LARGE  TO  HISTOGRAM') 

230  CONTINUE 
1  CONTINUE 

IF(NAS.LT.  2) STOP 
WRITE (6, 908) 

DO  330  IHIST-1, MAXSUB 
BOT-FLOATC IHIST-1 ) * . 01 
TOP-FLOAT ( IHIST) * . 01 
NPIII ST-NOTOTP  (IHIST) 

IFCNPHIST  .LE.  110) GO  TO  325 
WRITE (6, 910)  30T ,  TOP 
GO  TO  330 
325  CONTINUE 

IFCNPHIST. GT.O)WRITE(6,909)BOT, TOP, NPHIST, (IASTER, I AST-1, 
1NPHIST) 

909  FORMAT  ( IX , F5 . 2  ,F7 . 2 ,15 , 2X  ,1111 , 110A1 ) 

910  F0RMAT(1X,F5.2,F7.2,7X,'I  TOO  LARGE  TO  HISTOGRAM') 

330  CONTINUE 

STOP 

901  FORMAT (6F12. 5) 

902  FORMAT (613) 

904  FORMAT (3F10. 4) 

905  FORMAT ( 1H1 , 4X , 1HI , 5X , 1HH , 5X , 1HK , 5X , 1HL , 5X , 1HN , 7X , 4HDIST) 

906  FORMAT (5(1X,15) ,2X,F10.4,1X,A1) 

END 


B  Theoretical  EXAPS  Program 


Q 

C  THEORETICAL  EXAPS  PROGRAM  USING  EQ(2-3^) 

C 

REAL  K,K1 
REAL  C(*D 

DIMENSION  CHI2 (  501,30)  ,  v  , 

DIMENSION  ETA(  501),CHI1(  501), E(  501). K( 
DIMENSION  AN (  30) ,D( 30 )  , Phi ( 30) ,RH0( 30 ) 
DIMENSION  AP(i»),BF(4) 

DIMENSION  HOL(20) 

PI-3. 1^159265 
CMIN-1.E70 
CMAX-1.E70 
C 

C  LATTICE  PARAMETER  POR  COPPER 

C 

A-3.61 

c 

C  NS  IS  THE  NUMBER  OP  SHELLS  TO  BE  SUMMED 
C 

MS-13 

C 

C  READ  INPUT  DATA 

C 

READ (5,1)  HOL 

1  FORMAT (20 A4) 

READ(5,8)  (AP(I)  ,BF(I)  ,1-M) 

8  FORMAT ( 8F9 .  ll ) 

READ( 5 ,10 )  (AN(l),D(I) , PHI (I), 1=1, NS) 

10  FORMAT (3Fb.O) 

WRITE  6, 2) HOL 

2  FORMAT (1H1,20AH) 

WRITE (6 ,11)  (AN(I) ,PHI (I) ,1=1, NS) 

11  FORMAT (IX, 3*6. H) 

N = 4  9  9 

C 

C  GENERATE  ENERGY  VALUES 

C 

DO  100  1=1, N 
XI=I 

E(I)=2.*XI 
E1=E (I ) 

C 

C  DEFINE  K  USING  EQ(2-l6)  AND(2-35) 

C 

K(I)=SQF.T(.263*(E1+10)) 

K1=K(I) 


501) 


SS-Kl#Kl/(l6#PI#PI ) 

C 

C  DOYLE  AND  TURNER  SCATTERING  FACTOR 

C 

F=AF(1)»EXP(-BF(1)#SS)+AF(2)#EXP(-BF(2)#SS)+AF(3)*EXP(-BF(3) 
1  +AF(4)*EXP(-BF(4)*SS) 
c 

c  DEFINE  LINEAR  PHASE  SHIFT 

C  . 

ETA(I)=-.2i*5*Kl+7.95 

CHI-0.0 

c 

c  SUM  OVER  COORDINATION  SHELLS 

C 

DO  200  ,T=1  ,NS 
R=A*D( J ) 

ARG=K1»A#D(J) 

ARO>AR(HETA(I) 

RHO( J ) =-Kl*AN ( J ) *SIN( 2 . * ARG2)/ARG**2 
C 

C  SCATTERING  CORRECTION 

C 

RHO( J)=RHO(J )#K1*F 
C 

C  TEMPERATURE  CORRECTION 

C 

RH0( J)-RHO( J ) #EXP (-. 0265*K1**2*PHI ( J) ) 

C 

C  RANGE  TERM 

C 

RH0(J)«RH0(J)»KXP(-.3«D(J)»A) 

CIII=CHI+RHO(J) 

CHI2  (I ,  J )=CIII 

CHIN- AMIN1 ( CHIN , CHI 2 ( I , J ) ) 

CMAX-AMAX1  (CMAX  ,CHI2  ( IM.T) ) 

C(l)-CMIN 
C(2)-CMAX 
200  CONTINUE 
C 

C  WRITE  OUTPUT  DATA 
C 


30 

100 

C 


WRITE  (6,30)  el  ,K1  ,CH.I2  (1 ,1) ,  CHI  2  (1, 3)  CHI2(I  NS ) 
FORMAT  ( IX  ,F7 . 1  ,F8  i  ^  ,  3  (2X,E12 .  ^  ’  UsNo) 

CONTINUE 


C  PLOTTING  DATA  FOR  CALCOMP 
C 


INTEGER  I  BUF(102^ ) 
WRITE (6, 3) 

3  FORMAT( 11HOBEGIN  PLOT) 
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CALL  PLOTS (IBUF, 1024,0) 

CALL  PL0T(0. ,1.5, -3) 

CALL  SCALE (E, 10. ,N,1) 

CALL  SCALE (C,9.,2,l) 

CALL  AXIS ( 0 . , 0 . , 1HE , -1 , 10 . , 0 . ,E (N+L) ,E ( N+2 ) ) 
CALL  AXIS(0.,0.,3HCHI,3,9.,90.,C93),C(H)) 
CALL  SYMBOL(0 . 8 , . 1 ,HOL ,0 .  ,80) 

DO  iJO  «T=NS  ,NS 
CHI2 (N+l , J ) BC ( 3 ) 

CHI2(N+2,J)=C(4) 

CALL  LINE (E, CHI2(1,J),N, 1,0,0) 

40  CONTINUE  * 


4 


CALL  PLOT(12. ,-1.5,999) 
WRITE (6, 4) 

FORMAT (9HOEND  PLOT) 


STOP 

END 
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